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' Abstract. This is a review of the Riemann-Hilbert approach to the large N asymptotics 

, in random matrix models and its applications. We discuss the following topics: random 

' matrix models and orthogonal polynomials, the Riemann-Hilbert approach to the large N 

(~| , asymptotics of orthogonal polynomials and its applications to the problem of universality in 

^ ' random matrix models, the double scaling limits, the large N asymptotics of the partition 

function, and random matrix models with external source. 

\o 

(N 

r^'. 1. Introduction 

Ph. 

This article is a review of the Riemann-Hilbert approach to random matrix models. It is 
based on a series of 5 lectures given by the author at the miniprogram on "Random Matrices 
and their Applications" at the Centre de recherches math ematiques, Universite de Montreal, 
in June 2005. The review contains 5 lectures: 

^ ! Lecture 1. Random matrix models and orthogonal polynomials. 

I Lecture 2. Large asymptotics of orthogonal polynomials. The Riemann-Hilbert 

approach. 

00 ■ Lecture 3. Double scaling limit in a random matrix model. 

Lecture 4. Large N asymptotics of the partition function of random matrix models. 
Lecture 5. Random matrix models with external source. 

o ■ 

(yQ \ The author would like to thank John Harnad for his invitation to give the series of lectures 
O ■ at the miniprogram. The lectures are based on the joint works of the author with several 
coauthors: Alexander Its, Arno Kuijlaars, Alexander Aptekarev, and Bertrand Eynard. The 
^ ! author is grateful to his coauthors for an enjoyable collaboration. 

b : 

Lecture 1. Random matrix models and orthogonal polynomials 

The first lecture gives an introduction to random matrix models and their relations to 
orthogonal polynomials. 

2. Unitary ensembles of random matrices 

2.1. Unitary ensemble with polynomial interaction. Let M = {Mj^)^^^^ be a random 
Hermitian matrix, Mkj = Mj^, with respect to the probability distribution 

d^^{M) = -Le-^Ti-^WrfM, (2.1) 
Zjsf 



Date: June 26, 2008. 

The author is supported in part by the National Science Foundation (NSF) Grant DMS-0354962. 

1 



2 PAVEL M. BLEHER 

where 

p 

V{M) = Y,tjM^, P = 2po, tp>0, (2.2) 

i=l 

is a polynomial, 

N N 

dM ^Yl dMjj Yl dRe Mjkdira Mjk, (2.3) 

j=i j^k 

the Lebesgue measure, and 

Zjv=/" e-^^^(^)dM, (2.4) 

JHn 

the partition function. The distribution ixi^{dM) is invariant with respect to any unitary 
conjugation, 

M U-^MU, U e U{N), (2.5) 
hence the name of the ensemble. 

2.2. Gaussian unitary ensemble. For V{M) = M^, the measure //jv is the probability 
distribution of the Gaussian unitary ensemble (GUE). In this case, 

TV N 

TtV{M)=TtM^= 5]Mfe,M,-, = 5]M|,. + 2j]|M,fe|^ (2.6) 

j,k=l j=l j>k 

hence 

/^^""(^^) = ^ n (^"''''^) n (e-^^l^-l^) dM, (2.7) 

^ 3=1 j>k 

so that the matrix elements in GUE are independent Gaussian random variables. The 
partition function of GUE is evaluated as 

^s- = / n (-™^) n (--'""'0 « - iT' (w)""""' 

i=i i>k (2.8) 



2^ 

If V{M) is not quadratic then the matrix elements Mjk are dependent. 

2.3. Topological large N expansion. Consider the free energy of the unitary ensemble 
of random matrices, 

- -N-HuZn = -Ar-2 /" e-^^^(^)dM, (2.9) 

and the normalized free energy, 

F. = -7V-Mn^ = -iV-^ln ^Jj^^,^^^^ . (2.10) 

The normalized free energy can be cxpessed as 

F^ = -iV-2ln(e-^^^^W> (2.11) 
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where Vi{M) = V{M) - and 

(fiM)) - \ , (2.12) 

the mathematical expectation of f{M) with respect to GUE. Suppose that 

V{M) = + t3M^ + ... + tpMP. (2.13) 

Then (CTHl reduces to 

Fj, = _iv-2in(^e-^T^(*^*^'+-+*^^^^)) . (2.14) 
Fjv can be expanded into the asymptotic series in negative powers of A^^, 

9=1 

which is called the topological large N expansion. The Feynman diagrams representing F^^^^ 
are realized on a two-dimensional Riemann closed manifold of genus g, and, therefore, F/v 
serves as a generating function for enumeration of graphs on Riemannian manifolds, see, 
e.g., the works [BIZ], [BIPZ], [IZ], [DiF], [EM], [EMP]. This in turn leads to a fascinating 
relation between the matrix integrals and the quantum gravity, see, e.g., the works [DGZ], 
[Wit], and others. 

2.4. Ensemble of eigenvalues. The Weyl integral formula implies, see, e.g., [Mehj . that 
the distribution of eigenvalues of M with respect to the ensemble fiN is given as 



N 



dfx^iX) = i n(^^- - ^^)' n e-^^^'^^c^A, (2.16) 
j>k j=i 

where 

N 

Zn= Y[{\j-\kYYle-^^^^^^d\, dX = d\i...dXN. (2.17) 

j>k j=l 



Respectively, for GUE, 



1 ^ 
^^'^''''W = - y^kf n ^"""'^'^A, (2.18) 

N jyk j=l 



where 

TV 



GUE 



j>k j 

The constant Z^^^ is a Selberg integral, and its exact value is 

(27r)^/2 ^ 



/ I[(^j~h?l[e-'''UX. (2.19) 

i>k .7=1 
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see, e.g., |Meh] . The partition functions and are related as follows: 

^^-^N _ - Y[n\ (2.21) 



7 7GUE 1 ^ 



n=l 



One of the main problems is to evaluate the large N asymptotics of the partition function 
Zn and of the correlations between eigenvalues. 
The m-point correlation function is given as 

RmN{Xi, . . . ,Xm) = Jirz rr / Piy{Xi, . . . , Xiy)dXm+l ■ ■ ■ dxjy, (2.22) 

{N - J^N-m. 

where 

N 

Pn{xu ...,xm) = ZJ,^ \{{x, - Xkf n e-^^(^^). (2.23) 

j>k j=l 

The Dyson determinantal formula for correlation functions, see, e.g., |Meh] . is 

RmN{Xl, ...,Xm)= dct {Kj^^Xk, Xi))2i=i , (2.24) 

where 

N-l 

KN{x,y) = J2Mx)My) (2.25) 



n=0 

and 

M^) = ^Pn(x)e-^^(^)/^ (2.26) 
where Pn{x) = + a„_ix"~^ + . . . are monic orthogonal polynomials, 

P„(x)P^(x)e-^^(^')rfx = hJnm. (2.27) 

Observe that the functions ipn{x), n = 0, 1, 2, . . ., form an orthonormal basis in L^(M^), and 
Kjsi is the kernel of the projection operator on tne dimensional space generated by the 
first functions ipni n = 0, . . . , — 1. The kernel Kjsi can be expressed in terms of ipN-i, 
ipN only, due to the Christoffel-Darboux formula. Consider first recurrent and differential 
equations for the functions ijjn- 

2.5. Recurrence equations and discrete string equations for orthogonal polyno- 
mials. The orthogonal polynomials satisfy the three term recurrent relation, see, e.g. [Sze], 



In = 7 > 0, n > 1; 70 = 0. 



xPn{x) = Pn+lix) + (3nPn{x) + ^lPn-l{x), 
hn-1 / 

For the functions ipn it reads as 

X1pn{x) = 'Jn+lAi+lix) + Pn^n{x) + Jntpn-lix). (2.29) 
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This allows the following calculation: 



N-l 



N-1 



(X 



n=0 



n=0 



-ll^nix) {'Jn+l^n+l{y) + Pn^niv) + 7nV'n-l (Z/) )] 

= In [ipN{x)iJN~i{y) - i'N-i{x)i'Niy)] 



(2.30) 



(telescopic sum), hence 



N-l 



KN{x,y) = ^^„(x)^„(?/) 



7n- 



ipN{x)i)N-i{y) - i'N-i{x)i'Niy) 



n=0 



x-y 



which is the Christoffel-Darhoux formula. For the density function we obtain that 

Qn{x,x) In V n 



Pn[x) 



[iIj'^{x)^n~i{x) - ij'N^i{x)ijN{x)] 



(2.31) 



(2.32) 



N 

Consider a matrix Q of the operator of multiplication by x, f{x) — > xf{x), in the basis 
{ipnix)}- Then by fl2.29p . Q is the symmetric tridiagonal Jacobi matrix, 

//3o 71 ...\ 
7i A 72 ... 
Q= 72 /52 73 . . . . (2.33) 

73 /?3 ••• 

V: ; ; ; 

Let P = {Pnm)nm=oi2 ^ matrix of the operator f{z) — * f'{z) in the basis ipn{.z), 
n = 0,1,2, ... , so that 

r-OO 



Then P„ 
hence 



'Pnm. and 



NV'iz] 



Pr, 



PLiz 



i^n{x)i)'^{x) dx. 



(2.34) 



n\fJ_^~NV{z)/2 
hr. 



NV'iz 



n 



^„(^) + _V;„_i(z) + ..., (2.35) 

2 7n 



p + 
p + 



NV'{Q) 
2 

w(g) 



nn 

n 



P + 



w(g) 



0, 



n,n+l 



Since P^n = 0, we obtain that 



[V\Q)]nn = 0. 



In addition, 



hence 



-P + 



W(Q) 



n.n—l 



P + 



NV'iQ) 



n 

n,n— 1 7n 



7n[^'(Q)]n, 



n-l 



n 
iV 



(2.36) 

(2.37) 
(2.38) 
(2.39) 



n=l 



Then equations i2.4C^ can he written 



as 



71 = ■ (2.47) 



e 

oo 



The Hamiltonian is 



n=l 



2 7^ (2t + gil^i + gil + gil+i) -n\n^l 



(2.40) 
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Thus, we have the discrete string equations for the recurrent coefficients, 

7n[nQ)kn-l = ^, 
[V'mnn = 0. 

The string equations can be brought to a variational form. 
Proposition 2.1. Define the infinite Hamiltonian, 

oo 

H{',,D) = NTiV{Q)-Y,'^\n',l, 7 = (7o,7i, • • 0. /3 = (A, ft. . . .)■ (2.41) 



1^ = 0, 11^ = 0; „>1, (2.42) 

which are the Euler- Lagrange equations for the Hamiltonian H . 
Proof. We have that 

iVTr ( V'{Q)^\ -- = 2N[V'{Q)U^^, - -, (2.43) 



^'~1n V C'7n / 7rt 7n 

and 

hence equations fl2.40p are equivalent to (12.421) . □ 
Example. The even quartic model, 

V{M) = ^M^ + (2.45) 

In this case, since V is even, /?„ = 0, and we have one string equation, 

7n(^ + ^?7' 

with the initial conditions: 7o = and 



il [t + ^?7Li + gil + gil+i) = ^ , (2.46) 



(2.48) 



The minimization of the functional if is a useful procedure for a numerical solution of the 
string equations, see [BDJT], [BI2]. The problem with the initial value problem for the string 
equations, with the initial values 70 = and (I2.47p . is that it is very unstable, while the 
minimization of H with 70 = and some boundary conditions at n = A^, say = 0, works 
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very well. In fact, the boundary condition at n = creates a narrow boundary layer near 
n = N, and it does not affect significally the main part of the graph of 7,^. FigH] presents 
a computer solution, y = 7^, of the string equation for the quartic model: g = 1, t = — 1, 
= 400. For this solution, as shown in [BIl], [BI2], there is a critical value, Ac = so that 
for any e > 0, as N ^ 00, 



7^ = i?(^-j+0(iV-^), if ^>A, + e, (2.49) 
and 

[i?(^)+0(iV-^), n = 2k + l, 

[L(-)+0(iV-i), n = 2k, ^ 
The functions R for A > Ac and R, L for A < Ac can be found from string equation (12.461) : 



R{X) = ' + % + A>Ac, (2.51) 

and 

i?(A),L(A) = ii^^^Hi5, A<Ac. (2.52) 

We will discuss below how to justify asymptotics (12.491) . (12.501) . and their extension for a 
general V. 

2.6. Differential equations for the ^/'-f unctions. Define 
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Then 
where 



AJz) 



^'^{z) = iVA„(z)C(^), (2.54) 

(2.55) 



2 

y -'JnVn-l{z) — VlnUn[z)j 



and 
where 



Un{z) = [W{Q, z)]„,„_i, Vn{z) = [W {Q , z)]„„, (2.56) 

Q - z 

Observe that An{z) = 0. 

Example. Even quartic model, V{M) = fM^ + . Matrix An{z): 

"^^"V -7ni9z' + e^.,) \{tz + gz^)+g^lz) ^^'^^^ 

where 

en = t + g^l + g^l^,. (2.59) 

2.7. Lax pair for the discrete string equations. Three term recurrence relation (12.281) 
can be written as 

^n+i{z) = f/„(z)^„(z), (2.60) 

where 

Un{z) = (^''+^^\~ ^""^ -7n+l7n^ (2.61) 

Differential equation fl2.54p and recurrence equation fl2.60p form the Lax pair for discrete 
string equations fl2.40p . This means that the compatibility condition of (I2.54p and (12.601) . 

U'^ = N{An+lUn-UnAn), (2.62) 

when written for the matrix elements, implies fl2.40p . 

3. The Riemann-Hilbert problem for orthogonal polynomials 
3.1. Adjoint functions. Introduce the adjoint functions to Pn{z) as 

, 2GC\M, (3.1) 



°° PJu)w(u)du 



where 

w{z) = e-^^(^) (3.2) 
is the weight for the orthogonal polynomials P„. Define 

Qn±{x) = lim Qn{z), — OO < X < oo. (3.3) 

±lm^>0 

Then the well-known formula for the jump of the Cauchy type integral gives that 

Qn+{X) - Qn-{X) = w{x)Pn{x). (3.4) 
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The asymptotics of Qn{z) as z ^ oo, z E C, is given as 



^ , , I r w(u)Pju) du 

Qn[z) = — 



27rz u — z 



, oo 



(3.5) 



j=n+2 



(due to the orthogonahty, the first n terms cancel out). The sign = in (13.51) means an 
asymptotic expansion, so that for any k > n + 2, there exists a constant Ck > such that 
for all z E C, 

It can be some doubts in the uniformity of this asymptotics near the real axis, but we assume 
that the weight w{z) is analytic in a strip {z : |Im2;| < a}, a > 0, hence the contour of 
integration in (13.51) can be shifted, and (13. 6p holds uniformly in the complex plane. 

3.2. The Riemann-Hilbert problem. Introduce now the matrix-valued function, 

V ( ^\ — I Pn{z) Qn{z) \ , _^ 

^ " V^Pn-l(^) CQn^riz))^ ^^-'^ 

where the constant, 

C=-^, (3.8) 

is chosen in such a way that 

CQ„_i(z) = ^ + ..., (3.9) 

see (13.51) . The function y„ solves the following Riemann-Hilbert problem (RHP): 

(1) Yn{z) is analytic on C+ = {Im^ > 0} and C~ = {Im^ < 0} (two-valued on R = 

c+nc-). 

(2) For any real x, 

y;+(x) = F„_(x)jy(x), jy(x) = (^J (3.10) 

(3) As 2; ^ oo, 

^"(--)^(^+|:5)(o A) (3.11) 

where Y^, = 1, 2, . . . , are some constant 2x2 matrices. 

Observe that fl3J0D follows from ([33!), while fl3lTl) from ([33]). The RHP (l)-(3) has some 
nice properties. 

First of all, (13.71) is the only solution of the RHP. Let us sketch a proof of the uniqueness. 
It follows from fl3A0D . that 

det Yn+{x) = det F„_(a;), (3.12) 
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hence detYn{z) has no jump at the real axis, and hence Yn{z) is an entire function. At 
infinity, by (13.111) . 

detYniz) = l + ... (3.13) 

hence 

detYn{z) = l, (3.14) 

by the Liouville theorem. In particular, Yn{z) is invertible for any z. Suppose that Yn is 
another solution of the RHP. Then X„ = YnY~^ satisfies 

X„+(x) = Y„.4x)Yn+{x)-' = Yn4x)jy{x)jy{x)-'Y-_}{x) = X„_(a;), (3.15) 

hence X„ is an entire matrix- valued function. At infinity, by (13. lip . 

Xn{z)=I + ... (3.16) 

hence X„(z) = /, by the Liouville theorem. This implies that Yn = Yn, the uniqueness. 
The recurrent coefficients for the orthogonal polynomials can be found as 

7^ = [Fi]2i[ri]i2, (3.17) 



and 



f^n-i = - [Y^]u . (3.18) 



Indeed, from (ETD, (EUl), 

rn[Z) 1^ Q ^nj- \Cz-^Pn.liz) Z^CQn-l{z) i " ' + 2^ ' ^^'^^^ 

fc=l 



hence by (ES]), (ES}, and (ESHD, 



(3.0) 



which proves (13.171) . Also, 

[Y2U 



\Yl]ll = Pn-l,n-2 - Pn,n-1 , (3.21) 



where 

n 

Pn{z) = Y,PnjZK (3.22) 

i=o 

From (I2.28P we obtain that 

Pn-l,n-2 — Pn,n-1 = Pn-1 , (3.23) 

hence (13.180 follows. The normalizing constant /i„ can be found as 

27ri 

hn = -2n^[Y^]^2, K^^ = -—— . (3.24) 

The reproducing kernel K]s[{x,y) of the eigenvalue correlation functions, see (12.241) . is 
expressed in terms of Ym+{x) as follows: 

i^;v(a:,y) = e-^e-^^-^^— y(0 l) F4(y)F^+(x) Q . (3.25) 
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Indeed, by ( ICTjl . 

i^N{x)iJN-i{y) - i'N-i{x)iJN{y) 



(3.26) 



KN{x,y) = 7Ar- 

x-y 

NVM NVjy) PN{x)PN-i{y) - PN-i{x)PN{y) 

= e 2 e 2 = . 

\jhf<[hi^_i x — y 

From dXZD, dSSD and fIXTlD . we obtain that 

(0 1) Y^l{y)Y^^{x) L j = ^ [P^(a;)P^_i(y) - P^_i(a;)P^(y)] . (3.27) 



Also, 

7n 



(3.28) 



hence equation (I3.25P foUows. 

4. Distribution of eigenvalues and equilibrium measure 

4.1. Heuristics. We begin with some heuristic considerations to explain why we expect 

that the limiting distribution of eigenvalues solves a variational problem. Let us rewrite 
(EUD as 

d^M{X) = Z^^e~^^^^^dX, (4.1) 

where 

N 

Hn{X) = -J2^n\X,-Xk\+Nj2 Vi^j)- (4-2) 

j=i 

Given A, introduce the probability measure on M\ 

N 

dux{x) = N-^ 6{x - \j)dx. (4.3) 

i=i 

Then fl4.2p can be rewritten as 



Let u be an arbitrary probability measure on M}. Set 
/y(^) = - 

Then (14.1 p reads 



In \x — y\di'\{x)di>\{y) + / V{x)dv\{x) 



(4.4) 



// \n.\x — y\di'{x)di'{y) + I V{x)di'{x). (4.5) 

J J x^y J 



dfi^iX) = Z^^e-^''^('^^)rfA. (4.6) 

Because of the factor N"^ in the exponent, we expect that for large N the measure /xat is 
concentrated near the minimum of the functional Jy, i.e. near the equilibrium measure z/y. 
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4.2. Equilibrium measure. Consider the minimization problem 

Ev= inf Jy(z/), (4.7) 

!^eAfi(R) 

where 

Mi(M) = |z/ : z/ > 0, y di^ = l| , (4.8) 
the set of probability measures on the line. 

Proposition 4.1. (See [DKM].) The infinum of Iv{^) is attained uniquely at a measure 
V = vy , which is called an equilibrium measure. The measure vy is absolutely continuous, 
and it is supported by a finite union of intervals, J = U^^]^[aj, bj]. On the support, its density 
has the form 

Pv{x) = ~T—{x) = 7^ — :h{x)R^_/'^{x), R{x) = Y\i^ ^ aj){x — bj). (4.9) 



dx 27ii 



1/2, 



X] 



Here i?^/^(x) is the branch with cuts on J, which is positive for large positive x, and R 
is the value of R^^'^{x) on the upper part of the cut. The function h{x) is a polynomial, which 
is the polynomial part of the function O't infinity, i. e. 

= hix) + 0{x-^ (4.10) 

In particular, deg h = deg V — 1 — q. 

There is a useful formula for the equilibrium density [DKM]: 



(iz/y(x) 



-V^, (4.11) 



dx 71 

where 

This, in fact, is an equation on q, since the right-hand side contains an integration with 
respect to uy- Nevertheless, if \^ is a polynomial of degree p = 2pQ, then (14.121) determines 
uniquely more then a half of the coefficients of the polynomial q, 

q{x) = -(^^y -Oix^-'). (4.13) 

Example. If V{x) is convex then uy is regular, and the support of uy consists of a 
single interval, see e.g. |KuMlj . For the Gaussian ensemble, V{x) = x^, hence, by f l4.13p . 
q{x) = — x"^. Since 

•-^ 1 

— va^^—x^dx = 1, 

we find that a = a/2, hence 

py{x) = -^2 - x2 , \x\<V2, (4.14) 

TT 

the Wigner semicircle law. 
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4.3. The Euler-Lagrange variational conditions. A nice and important property of 
minimization problem (14. 7p is that the minimizer is uniquely determined by the Euler- 
Lagrange variational conditions: for some real constant /, 

2 / log \x - y\diy{y) - V{x) = I, for x E J, (4-15) 
2 / log \x - y\dv{y) - V{x) < /, for x G M \ J, (4.16) 

see [DKM]. 

Definition. (See [DKMVZ2].) The equilibrium measure, 

dvv{x) = ^h{x)R^''^{x)dx (4.17) 
27rz 

is ca/Zed regular ('ot/ierazse singular^ if 

(1) h{x) ^ on the (closed) set J. 

(2) Inequality (^T7^ is strict, 

2 j log \x - y\dvv{y) - V{x) <l, for X G M \ J. (4.18) 

4.4. Construction of the equilibrium measure: equations on the end-points. The 

strategy to construct the equilibrium measure is the following: first we find the end-points 
of the support, and then we use equation (14.101) to find h{x) and hence the density. The 
number q of cuts is not, in general, known, and we try different g's. Consider the resolvent, 

u{z)= .GC\J. (4.19) 

Jj z-x 

The Euler-Lagrange variational condition implies that 



.W = m_M£)«!^!(£). ,4.20) 



Observe that as 2; ^ cxd, 



u;{z) = - + ^ + ..., mk= x^p{x)dx. (4.21) 

z Z J J 



The equation 

v'(.) ^^^^^^ 



2 2 z 

gives g + 1 equation on oi, 61, . . . , a^, 6^, if we substitute formula (]4.10p for h. Remaining 
g — 1 equation are 

'^\{x)R^/'^{x)dx = Q, j = l,...,g-l, (4.23) 



which follow from (I4.20p and (14.151) . 

Example. Even quartic model, V{M) = + \M^. For t > tc = —2, the support of 
the equilibrium distribution consists of one interval [—a, a] where 



l/9\ 1/2 

-2t + 2(t2 + l2)'/'^ 



(4.24) 
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Figure 2. The density function, pv{x), for the even quartic potential, 
V{M) = + iM^ for t = -1, -2, -3. 



and 



where 



Pv{x) 



vr V 2 ^ 

t + ((tV4) + 3)'/' 



In particular, for t = —2, 



Pvix) = — a;^V4 — a;^ 
27r 



For t < —2, the support consists of two intervals, [—a, —6] and [6, a], where 

a = V2 - t, 6= V-2-t, 
and ^ 

— |x| (a^ — x'^){x'^ — 6^) . 
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(4.25) 

(4.26) 

(4.27) 

(4.28) 
(4.29) 



Figl2] shows the density function for the even quartic potential, for t = —1, —2, —3. 



Lecture 2. Large N asymptotics of orthogonal polynomials. The 
Riemann-Hilbert approach 

In this lecture we present the Riemann-Hilbert approach to the large asymptotics of 
orthogonal polynomials. The central point of this approach is a construction of an asymptotic 
solution to the RHP, as — ^ oo. We call such a solution, a parametrix. In the original paper 
of Bleher and Its [BIl] the RH approach was developed for an even quartic polynomial V{M) 
via a semiclassical solution of the differential equation for orthogonal polynomials. Then, in a 
series of papers, Deift, Kriecherbauer, McLaughlin, Venakides, and Zhou [DKM], [DKMVZl], 
[DKMVZ2] developed the RH approach for a general real analytic V, with some conditions 
on the growth at infinity. The DKMVZ-approach is based on the Deift-Zhou steepest descent 
method, see |DZ]. In this lecture we present the main steps of the DKMVZ-approach. For 
the sake of simplicity, we will assume that V is regular. In this approach a sequence of 
transformations of the RHP is constructed, which reduces the RHP to a simple RHP which 
can be solved by a series of perturbation theory. This sequence of transformations gives the 
parametrix of the RHP in different regions on complex plane. The motivation for the first 
transformation comes from the Heine formula for orthogonal polynomials. 
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5. Heine's formula for orthogonal polynomials. 

The Heine formula, see, e.g., [Sze], gives the N-th orthogonal polynomial as the 
integral, 



Pn(z) 



{det{z -M)) = Z^'I' det{z - M)e~^^'^'^^'^dM. 



In the ensemble of eigenvalues, 

N \ „ N N 

pn(z) = ^ - ^^•) i^^N - ^^•) n(^^- - ^^)' n ^-^^^^'^^^a 

7 = 1 / 1 = 1 j>k 1=1 



TV . 

\i=i 



Since ux is close to the equilibrium measure u for typical A, we expect that 



\j=i 

hence by the Heine formula. 



N~' log (y[{z - A,)\ ^ \og{z - x)dMx), 



N ^\ogPN{z)^ J \og{z — x)diyv{x). 
This gives a heuristic semiclassical approximation for the orthogonal polynomial. 



/ log(2 — x)di'v{x) 



Pn{z) ^ exp 

and it motivates the introduction of the "gf-function" . 
5.1. (yf-function. Define the (^-function as 

g{z) = J \og{z — x)diyvix), zeC\{—oo,bq] 

where we take the principal branch for logarithm. 
Properties of g{z): 

(1) g{z) is analytic in C \ (— oo, bg]. 

(2) As z ^ oo 

g{z) = \ogz -y2^ ^ 9j= / —duv{x). 

(3) By dnni), g^OD, 

(4) By dUSl), 

g+{x) + g^{x) = V{x) + /, x e J. 

(5) By iH]), 

g+{x) + g^{x) <V{x) + 1, xeR\J. 



= 27Tpv{x) > 0, X e {ttj, bj), 1 < j <q. (5-14) 
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(6) Equation (15.51) implies that the function 

G{x) = g^x) - g4x) (5.10) 
is pure imaginary for all real x, and G{x) is constant in each component of M \ J, 

G{x) = iVtj for bj < X < a^+i, 1 < j < q — 1, (5-11) 

where 

= 27r ^ / pvix) dx, 1 < j < g - 1. (5.12) 
k=j+i "^"fc 

(7) Also, 

/•X 

G{x) = iVLj — 2m / pv{s) ds for aj < x < bj, 1 < j < (5.13) 

where we set = 0. 
Observe that from (15.131) and (14.91) we obtain that G{x) is analytic on {aj,bj), and 

dG{x + iy) 
dy 

From (15. 8p we have also that 

G{x) = 2g+{x)-V{x)~l = -[2g_{x)-V{x)-l], x e J. (5.15) 

6. First Transformation of the RH Problem 

Our goal is to construct an asymptotic solution to RHP (I3.10p . (13. lip for Y^i^z), as N ^ 
oo. In our construction we will assume that the equilibrium measure is regular. By (15.40 
we expect that 

P^(z)^e^^(^), (6.1) 
therefore, we make the following substitution in the RHP: 

F^(z)=ef-T^(.)eM^(^)-^]-, ^3 = . (6.2) 

Then Tn{z) solves the following RH problem: 

(1) Tn{z) is analytic in C \ M. 

(2) Tn+{x) = T]s[-{x)jT{x) for x G M, where 

/^-N[g+{x)-9.{x)] ^N[9+{x)+g^{x)-V(x)-l]\ 
jrix) = i Q ^N[g+ix)-g.{x)] J ■ (6-3) 

(3) Tn{z) = I + 0{z-^), as 2 ^ oo. 

The above properties of g{z) ensure the following properties of the jump matrix Jt'- 

(1) jrix) is exponentially close to the identity matrix on (—00, ai) U (6^, 00). Namely, 

Jr(a;) = L M, x G (-00, ai) U (6^, 00), (6.4) 
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where c{x) > is a continuous function such that 

c{ X] 

hm - — ■ — - = oo, hm c{x) = hm c{x) = 0. (6-5) 

X— utoo In X x^ai x-^bq 



(2) For 1 < J < g - 1, 



iT{.x)=y Q ^^iNn, y xe{bj,aj+i), (6.6) 

where c{x) > is a continuous function such that 

hm c{x) = hm c(x) = 0. (6.7) 

x—fbj a;— >aj+i 

(3) On J, 

3t{x) = ( Q g7VG(x) 1 • (6-8) 

The latter matrix can be factorized as foUows: 

'g-JVG(x) ^ \ f I 0\/Ol\/l 



^NG{x) j - y^NG{x) oy^e~^G(-) l) (6.9) 

This leads to the second transformation of the RHP. 

7. Second transformation of the RHP: Opening of lenses 

The function Q-^^^i^) is analytic on each open interval (a-,-, hj). Observe that [e"^*^*^^)! = 1 
for real x E {aj,bj), and e"^*^^^^ is exponentially decaying for Imz > 0. More precisely, by 
f l5.14p . there exists yo > such that e"^'^^^) satisfies the estimate, 

|g-7VG(.)| < e-Nc{z)^ ^^^^^^^y. aj <x <bj, 0<y < yo}, (7.1) 

where c{z) > is a continuous function in R^. Observe that c{z) — > as Imz — > 0. In 
addition, |e^'^(^)| = |e~^*^(^^|, hence 

|gjvG(.)| < e'Nciz)^ z e RJ = {z = x + iy : aj < x < bj, < -y < yo}, (7.2) 
where c{z) = c{z) > 0. Consider a C°° curve 7^ from aj to bj such that 

1^ = {x + iy: y = fj{x)}, (7.3) 
where fj{x) is a C°° function on such that 

f,{aj) = f,{b,) = 0; /;(a,) = -/'(6,) = v^; < /,(x) < yo, aj < x < bj. (7.4) 
Consider the conjugate curve, 

77 =1^ = {x-iy : y = fj{x)}, (7.5) 

see FigJ3J The region bounded by the interval [aj, bj] and 7^ (7^) is called the upper (lower) 
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Figure 3. The lenses. 

lens, C^, respectively. Define for j = 1, . . . , g, 

TN{z)j^^{z), if z is in the upper lens, z E , 
Sn{z) — \ Tj^{z)j^{z), if z is in the lower lens, z e (7.6) 
Tn{z) otherwise. 



where 



3±{z) = ( ^:^NG{z) \ ] ■ (7.7) 



Then Sn{z) solves the following RH problem: 

(1) Sn{z) is analytic in C \ (R U T), T = 7^+ U 7^ U • • • U 7+ U 7". 

(2) 

Sr,+iz)^SN4z)js{z), zeRU^, (7.8) 
where the jump matrix js{z) has the following properties: 

for z G J. 



(a) Jsiz) 




(b) js{z) 


= jriz) 




= jriz) 


(c) js{z) 


= J±{z) 



-iNU 



Q \iNn, j for ^ e [bj, aj+i), j = 1, 1, and 
(q M for ^ e (-00, ai) U (65,00). 

(^Q(^^^c{z)N^ ior z e -ff, j = 1,. .. ,q, where c{z) > is a 

continuous function such that c{z) — >• as ^ — >• a^, bj. 
(3) Sn{z) = I + 0{z~^), as z ^ 00. 

We expect, and this will be justified later, that as N ^ 00, Sn{z) converges to a solution 
of the model RHP, in which we drop the 0(e~'^^)-terms in the jump matrix jsiz)- Let us 
consider the model RHP. 



8. Model RHP 

We arc looking for M{z) that solves the following model RHP: 
(1) M{z) is analytic in C \ [ai, bg], 
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(2) 

M+{z) = M_{z)jm{z), z e [a^bg], (8.1) 
where the jump matrix juiz) is given by the following formulas: 

(a) jmiz) = for z e J 

(b) iM(^) = ( Q ^jjvnj for ^ e [6j,aj+i], j = 1,. . 1. 

(3) M{z) = 1 + 0{z-^), as z^oo. 
We will construct a solution to the model RHP by following the work [DKMVZ2] . 

8.1. Solution of the model RHP. One-cut case. Assume that J consist of a single 
interval [a, b]. Then the model RH problem reduces to the following: 

(1) M{z) is analytic in C \ [a, b], 

(2) M+iz) = M4z) for z e [a,b]. 

(3) M{z) = 7 + 0{z-^), as z^oo. 

This RHP can be reduced to a pair of scalar RH problems. We have that 



A _ 1 A l \ fi \ fl -i 
-1 Oj ~ 2 [i -i [O -i \l i 



(2) M+(^) = M_(^)(^^ for [a, 6]. 

(3) M{z) = 7 + 0{z-^), asz^oo. 



I fl l \ f^-^ 0\ A 1 



detM(z) = 1. 
At infinity we have 



(8.2) 



Let 

MW = i(; -f)Mi^(\ (8.3) 

Then M{z) solves the following RHP: 
1^1) M{z) is analytic in C \ [a, b], 



This is a pair of scalar RH problems, which can be solved by the Cauchy integral: 
where 

7(.) - (^) (8.5) 
with cut on [a, b] and the branch such that 7(00) = 1. Thus, 



^^^^ 2 V -ij V Tyl^i -ij y iiz)-i-\z) lizj+j-^z) J ' (g g) 



7(^) = 1 + ^ + 0(^-2), (8.7) 
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hence 



/ n b-a 



z \^ 



MW=/ + -U„, ^-''n+0{z-'). (8.8) 



8.2. Solution of the model RHP. Multicut case. This will be done in three steps. 
Step 1. Consider the auxiliary RHP, 

(1) Q{z) is analytic in C \ J, J = Vfj^-^[aj, hj], 

(2) g+(^) = Q-(^)(_°i J) ioizeJ. 

(3) Q{z) = 1 + 0{z-^), as 2 ^ oo. 

Then, similar to the one-cut case, this RHP is reduced to two scalar RHPs, and the solution 
is 



where 



7(z)-7 ^(2) 7(z)+7 ^jz) 
2i 2 



. 1/4 



with cuts on J. At infinity we have 



=11(735^) ' 7(00) = 1, (8.10) 



7(z) = l + X:^ + 0(^-^), (8.11) 



Az 



hence 



1/0 



Q(^) = ^ + - ..-a, ^-''^ + 0(.-^). (8.12) 

In what follows, we will modify this solution to satisfy part (b) in jump matrix in (18. ip . This 
requires some Riemannian geometry and the theta function. 

Step 2. Let X be the two-sheeted Riemannian surface of the genus 

g = q-l, (8.13) 

associated to ^/R{z), where 



R{z) = l[{z-a,){z-bj), (8.14) 

j=i 

with cuts on the intervals {aj,bj), j = 1, . . . ,q, see FigJH We fix the first sheet of X by the 
condition that on this sheet, 

^/R(x)>0, x>bg. (8.15) 

We would like to introduce 2g cycles on X, forming a homology basis. To that end, consider, 
for j = 1, . . . , g, a cycle Aj on X, which goes around the interval {bj, Oj+i) in the negative 
direction, such that the part of Aj in the upper half-plane, A'j' ^ Aj U {z : Im^; > 0}, lies 
on the first sheet of X, and the one in the lower half-plane, Aj = Aj U {z : Imz < 0}, lies 
on the second sheet, j = 1, . . . , g. In addition, consider a cycle Bj on X, which goes around 
the interval {ai,bj) on the first sheet in the negative direction, see FigjSl Then the cycles 
{Ai, . . . , Ag, Bi, . . . , Bg) form a canonical homology basis for X. 
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X 



aj bj a 




Figure 4. The Riemannian surface associated to \fR{ 



z . 




I T", 



Figure 5. The basis of cycles on X. 



Consider the hnear space VL of holomorphic one- forms on X, 

The dimension of Vt is equal to g. Consider the basis in f2, 
with normalization 

/ uJk = Sjk, j,k = l,...,g. 

Ja, 

Such a basis exists and it is unique, see |FKj . Observe that the numbers 



z 



l<J<g, l<k<g-l, 



X) 



(8.16) 



(8.17) 
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are real. This implies that the basis uj is real, i.e., the one- forms. 



9-^ fc J 

CjkZ az 



3 



have real coefficients Cjk- 

Define the associated Riemann matrix of -B-periods 



(8-19) 



T = i^jk), Tjk = ujk, j,k = l,...,g. (8.20) 



Since ^JR{x) is pure imaginary on (a^, the numbers Tjk are pure imaginary. It is known, 
see, e.g., |FK] . that the matrix r is symmetric and {—ir) is positive definite. 
The Riemann theta function with the matrix r is defined as 

9 

meZf i=i 

The quadratic form i{m,Tm) is negative definite, hence the series is absolutely convergent 
and 6{s) is analytic in C^. The theta function is an even function, 

e{-s) = e{s), (8.22) 

and it has the following periodicity properties: 

e{s + ej) = e{sy, e{s + r,) = e-'^^''^-'''^^W{s), (8.23) 

where ej = (0, . . . , 0, 1, 0, . . . , 0) is the j-th basis vector in C^, and Tj = rcj. This implies 
that the function 

where c, c/ G are arbitrary constant vectors, has the periodicity properties, 

fis + e,) = fis); f{s + T,) = e-^^''^f{s). (8.25) 

Consider now the theta function associated with the Riemann surface X. It is defined as 
follows. Introduce the vector function, 

u{z)= [ iu, zeC\{ai,bg), (8.26) 

where u is the basis of holomorphic one- forms, determined by equations (18.171) . The contour 
of integration in (I8.26P lies in C \ (oi, bg), on the first sheet of X. We will consider u{z) as a 
function with values in C^/Z^, 

u : C\{ai,bg) CVZ^ (8.27) 

On [ai,6q] the function u{z) is two-valued. From (I8.20p we have that 

u+(x) — u^{x) = Tj, X G [bj, ctj+i]; 1 < j < — 1- (8.28) 

Since a/ R{x)_ = —\jR{x)_^ on [a^, bj], we have that the function u+{x) + U-{x) is constant 
on [aj,bj]. It follows from fl8.17p that mod 7J^, 

u+{bj) + u^{bj) = u+{aj+i) + u_{aj+i), 1 < j < g - 1 (8.29) 
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Since u+{hq) = U-{hq) = 0, we obtain that 

g 

u+{x) + u-{x) = 0, X e J =[J[aj,bj]. (8.30) 

j=i 

Define 

9{u{z) + d + c) 9{-u{z) +d + c) ^mx/ h\ (q 

fi{z) = / N , n > f2{z) = —^, . ^ , , zeC\{aubg), 8.31 

where c,d are arbitrary constant vectors. Then from fl8.28p and (18.251) we obtain that 
for 1 < j < g — 1, 

h{x + 20) = e-2--^/i(x - 20), h{x + zO) = e2--^/2(x - zO), x G (6„ a,+i), (8.32) 
and from (IHIBOj) that 

/i(x + iO) = /2(x-20), /2(x + i0) = /i(x-20), xeJ. (8.33) 

Let us take 

c = — , ^] = (^]l,...,^],), (8.34) 



27r 

and define the matrix-valued function 



(u(z)+di+c) e{-u{z)+di+c) \ 

_ I e(u(z)+di) e(-u{z)+di) ] /o OCX 

6l(«(^)+d2) 6l(-M(z)+d2) / 

where di,d2 G are arbitrary constant vectors. Then from (18.331) . we obtain that 
F+(x) = F_(x) ( Q ^i^f^, ) , X G {bj, aj+i); j = 1, . . . , g - 1, 



F^(x) = F_(x) ( J J 1 , xeJ. 



^ e{-u(oo)+d2) 



Then M(2;) has the following jumps: 



M+{x) = M_(x) ( Q ^ijvn, ) , a; e (6j, a^+i); j = 1, . . . , g - 1, 
M+(x) = M_(x)(_°^ XGJ, 



(8.36) 



Step 3. Let us combine formulae (18. 9p and (I8.35p . and let us set 

( 7(z)+7~Mz) 9{u{z)+di+c) 7(^)-7-i(^) 6>(-M(^)+di+c) \ 
2 e(u{z)+di) {-2i) e(~u{z)+di) ] /o oyN 

2i 6»(«(z)+d2) 2 e(-«(z)+d2) / 

where 

/ d{u(00)+dl+c) g \ 

i^(oo) = I '^-^T''^ ei-uioo)+d,+c) (8.38) 



(8.39) 



which fits perfectly to the model RHP, and M{oo) = I. It remains to find di, c?2 such that 
M{z) is analytic at the zeros of 6{±u{z) + di^2)- These zeros can be cancelled by the zeros 
of the functions 7(z) ± 7~^(z). Let us consider the latter zeros. 
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The zeros of ^{z) ± 7^^-^) are the ones of 7^(-2) ± 1, and hence of ■j^{z) — 1. By f lS.lOp . 
the equation 7"^ (2;) — 1 =0 reads 

P(^) = n ^ = 1 • (8.40) 

j=l J 

It is easy to see that 

p(6,+0) = oo, p(a,+i)=0, l<j<g-l, (8.41) 

hence equation fl8.40p has a solution xj on each interval (6j,aj+i), 

p(xj) = 1, bj < Xj < aj+i, 1 < j < g - 1. (8.42) 

Since equation fl8.4UI) has (g — 1) finite solutions, the numbers {xj, I < j < q — l} are all the 

solutions of fl8.40l) . The function 7(2;), defined by equation fl8.10p . with cuts on J, is positive 
on M \ J, hence 

7(x,) = 1. (8.43) 

Thus, we have (g — 1) zeros of 7(z) — j^^iz) and no zeros of 7(2;) + ^"^{z) on the sheet of 
7(2;) under consideration. 

Let us consider the zeros of the function 6{u{z) — d). The vector of Riemann constants is 
given by the formula 

q-l 

K = -J2<bj). (8.44) 

i=i 

Define 

q-l 

d= -K + (8-45) 
j=i 

Then 

9{u{xj) ~ d) = 0, l<j<q-l, (8.46) 

see [DKMVZ2], and {xj, l<j<g— 1} are all the zeros of the function 6{u{z) — d). In 

addition, the function 9{u{z) + d) has no zeros at all on the upper sheet of X. In fact, all 
the zeros of 6{u{z) + d) lie on the lower sheet, above the same points {xj, l<j<g— 1}. 
Therefore, we set in (18.370 . 

d^ = d, d2 = -d, (8.47) 

so that 

( 7(2)+7-i(z) e{u{z)+d+c) 7(z)-7-i(z) e{-u{z)+d+c) \ 

2 e{uiz)+d) {-2i) e{-u{z)+d) \ /o.ox 

7(z)--7-^(2) e{u{z)-d+c) ■y{z)+'y-^{z) e{-u{z)-d+c) I ' [O.^O) 
2i e{u{z)-d) 2 e{-u{z)-d) J 

where 

( 6>(M(00) + d+c) g \ 

'^"^^^^'^ .(-»(oo)-.+c) • (8.49) 

This gives the required solution of the model RHP. As 2; 00, 

M(z) = /+ — + 0(z^2)^ ^850) 
z 
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where 

Ml 







e(-u+d+c)e{u+d) \-^q (bj-aj) " 
e{u+d+c)9(-u+d) ^j=l (-4i) 
9{u—d+c)9(—u—d) ^ 
9\~u-d+c)e(u-d) 2^j=l ~ ^ 



U = UiOO 



1.51) 



Ai 



9. Construction of a parametrix at edge points 



We consider small disks D{aj, r), D{bj,r), 1 < j < of radius r > 0, centered at the edge 
points, 

D{a,r) = {z : l^: — a| < r}, 
and we look for a local parametrix U]\f{z), defined on the union of these disks, such that 
• Un{z) is analytic on D \ (R U F), where 

D = |J(D(a,,r)UZ}(6„r)). (9.1) 



Un+{z) = Un^{z)3s{z), ;2 g (m u r) n 



as — s> oo, 



Un{z) = [I + [ — ]] M{z) uniformly for z E dD. 



(9.2) 
(9.3) 




Figure 6. Partition of a neighborhood of the edge point. 



We consider here the edge point bj, I < j < q, in detail. From (l5.8p and (15.131) . we obtain 
that 



hence 



Pvis) ds, ttj < X < bj, 
[2g+{x) - V{x)] - [2g+{bj) - V{bj)] = -27ii [ pv{s) ds. 



(9.4) 



(9.5) 
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By using formula (14. 9p . we obtain that 

[2g^{hi)-V{h,)]-[2g+{x)-V{x)]= / h{s)Rf{s)ds, R{z) =\{{x - a,){x -h,). (9.6) 

Since both g{z) and R^^'^{z) are analytic in the upper half-plane, we can extend this equation 
to the upper half-plane, 

[2g4b,) - V{b^)] - [2g{z) - V{z)] = f h{s)R'/'{s) ds, (9.7) 

where the contour of integration lies in the upper half-plane. Observe that 

h{s)R^/^{s) ds = c{z - h,f/^ + 0{{z- hjf^) (9.8) 



as 2; — >■ 6j, where c > 0. Then it follows that 

m = {^[2^7+(&.) - V{h,)] - [2g{z) - V{z)]Y (9.9) 

is analytic at 6^, real-valued on the real axis near hq and I3'{hq) > 0. So /3 is a conformal 
map from D{bq,r) to a convex neighborhood of the origin, if r is sufficiently small (which 
we assume to be the case). We take F near bq such that 

P{T n D{bq, r)) C {z I arg(;2) = ±2n/3}. 

Then F and M divide the disk D{zi,r) into four regions numbered I, II, III, and IV, such 
that < arg/5(z) < 27r/3, 27r/3 < arg/3(z) < vr, -vr < aig l3{z) < -27r/3, and -27r/3 < 
SiTgP{z) < for z in regions I, II, III, and IV, respectively, see Figj6l 
Recall that the jumps js near bq are given as 

= J ) on [bj-r,bj) 

Js = I g-jvG(.) J on 7^ 

/ X (9.10) 

JS = L^G(.) J on 7g 

js=y ^N[g+{z)-g^(z)] j OU (6„ bj + v]. 

We look for Un{z) in the form, 

UNiz) = Q7vWe'^[^(^)-^-^>^ (9.11) 
Then the jump condition on Un{z), (19.21) . is transformed to the jump condition on Q^^Z), 

QN4z)=QN-iz)jQiz), (9.12) 

where 

jQiz) = e~^[^-(^)-^-^]-js(^)e^[^+(^)-'^-^]-. (9.13) 
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From (19A0D . (iBlSi) and (157[5|) we obtain that 

f 
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J) on - r, 6j) 



1 

1 1 

1 

1 1 

1 1 

1 



on 77 , 
on % , 

on {hj, bj + r]. 



(9.14) 



We construct Qn{z) with the help of the Airy function. The Airy function Ai{z) solves 
the equation y" = zy and for any e > 0, in the sector n + e < aigz < n — e, it has the 
asymptotics as 2; — > 00, 

1 



Ai(^) 



e-i^'^' (1 + 0(^-3/2)) _ 



(9.15) 



The functions Ai{ujz), Ai(to'^2;), where uj = e~ , also solve the equation y" = zy, and we 
have the linear relation, 

Ai(^) + ujAi{ujz) + uj^Ai{uj^z) = 0. (9.16) 

We write 

yoiz) = Ai{z), yi{z) = ujAi{ujz), y2{z) = uj^Ai{iu^z), (9.17) 
and we use these functions to define 

fyo{z) -2/2(2;)' 

yVoi^) -y'2{z), 

-y,{z) -y,{z)\ for27r/3<arg^<7r, 



, for < argz < 27r/3, 



$(z) = <^ 



-y'liz) -2/2 (^) 

-2/2(2) yi{z) 
-y'^{z) y[{z) 

'yo{z) yi{z) 

Observe that equation (19.161) reads 

2/0(2;) +2/1(2;) +2/2(2;) = 
and it implies that on the discontinuity rays, 



(9.18) 



for —TV < aigz < — 27r/3, 
, for — 27r/3 < arg^; < 0. 



2tt 



Now we set 
so that 



^+{z) = ^_{z)jq{z), arg2; = 0,±y ,7r. 



Qn{z)=En{z)<I>{N^/''P{z)), 
Un{z) = ^iv(2)$(iV'/'/5(^))e-'^[^(^)-^-5]'^^ 



(9.19) 

(9.20) 

(9.21) 
(9.22) 
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where is an analytic prefactor that takes care of the matching condition (19.31) . Since 
$(z) has the jumps jq, we obtain that Un{z) has the jumps js, so that it satisfies jump 
condition (19. 2p . The analytic prefactor i^^r is explicitly given by the formula, 

E^iz) = M{z)Qn{z)Ln{z)-\ (9.23) 

where M{z) is the solution of the model RHP, 

Q^[z) = e^ — ''\ ±lmz>0. (9.24) 

and 

1 [N-^l^ 13-^1 \z) W 1 



where for l3^/^{z) we take a branch which is positive for z G {bj,bj + r], with a cut on 
[bj — r, bj). To prove the analyticity of Ei<^{z), observe that 

[M{x)Qn{x)]+ = [M{x)Qn{x)]^Ji{x), bj - r <x <bj + r, (9.26) 

where 

j,(x) = e—'''3M{x)e—''\ (9.27) 

From (18. ip we obtain that 

iiix) = \ ^ z\ , bj - r < X <bj, 

^) ~ (9.28) 

ii{x) = I, bj < X < bj + r. 

From (K^ . 

L]sf+{x) = LN_{x)j2{x), bj-r<x<bj + r, (9.29) 
where j2{x) = I for bj < x < bj + r, and 

1 A f-i 0\ f 1 i\ / 1 



^■^W=(-l ^) [o = oj- (9-30) 

SO that j2{x) = ji{x), bj—r < x < bj+r. Therefore, E]\[[z) has no jump on bj—r < x < bj+r. 
Since the entries of both M and L have at most fourth-root singularities at bj, the function 
E]\f{z) has a removable singularity at z = bj, hence it is analytic in D{bj,r). 

Let us prove matching condition (19. 3p . Consider first z in domain I on FigEl From (l9.15l) 
we obtain that for < argz < 

^ . (9.31) 

hence for z in domain I, 

<l>(Ar2/3/5(2)) = ^ Ar-i-«/?(z)-3'^« ( \ (/ + 0(iV-i))e-i^^(")'^'"3 (9.32) 
From dlSD, 

^ = I {[294b,) - V{b,)] - [2g{z) - V{z)]} , (9.33) 
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• • • 




Figure 7. The contour Tr for Rn{z). 



q 



and from (]9.4p . 
hence 



2^+(&i) - = / + i^-, (9.34) 

- m'^' = ~g{z) + ^ + - + ^ . (9.35) 
Therefore, from (19.221) and (19.321) we obtain that 

1 1 1 / 1 A 1 

Um{z) = E^(^)_^iV-6-/3(^)-i- 1^^^ ^g_3g^ 

Then, from dOS]) and (lOSj) . 

?7^(^) = M(z)e^L^(^)-^-^ iV-6-3/?(^)-3-3 M (J + 0(^-1))^-^ 

2\/7r V (9.37) 

iiVQ ■ iiVQ ■ 

= M{z)e—{I + 0(iV-^))e-— = M(z)(/ + 0{N"^)), 

which proves (19. 3p for z in region I. Similar calculations can be done for regions II, III, and 
IV. 

10. Third and final transformation of the RHP 
In the third and final transformation we put 

Rn{z) = SNiz)M(z)~^ for z outside the disks D(aj,r), D{bj,r), I < j < q, 
Rn{z) = Sn{z)Un{z)~^ for z inside the disks. 

Then Rn{z) is analytic on C \ Tr, where Tr consists of the circles dD{aj,r), dD{bj,r), 
1 < J < Q', the parts of F outside of the disks D{aj,r), D{bj,r), 1 < j < q, and the real 
intervals (— oo, ai — r), {hi + r, 02 — r),. . . , (&g-i, ctg), {bg + r, 00), see FigJTl There are the 
jump relations, 

R^+{z) = RM-iz)jR{z), (10.2) 

where 

jniz) = M{z)Un{zY^ on the circles, oriented counterclockwise, 
jR{z) = M{z)js{z)M(z)~^ on the remaining parts of F^^. 



30 PAVEL M. BLEHER 

We have that 

Jr^z) = I + 0(N~^) uniformly on the circles, 

jniz) = I + 0{e~^^^^^) for some c{z) > 0, on the remaining parts of Tn. 
In addition, as x —>■ oo, we have estimate (I6.5p on c{x). As z —>■ oo, we have 



(10.4) 



°° R 

Rn{z) = I + J2^- (10.5) 

i=i 

Thus, Rn{z) solves the following RHP: 

(1) Rn{z) is analytic in C \ Tn. and it is two-valued on F/j. 

(2) On Tr, Rn{z) satisfies jump condition f ll0.2p . where the jump matrix jR{z) satisfies 
estimates f ll0.4p . 

(3) As z ^ cxD, Rn{z) has asymptotic expansion (110. 5p . 

This RHP can be solved by a perturbation theory series. 

11. Solution of the RHP for RNiz) 

Set 

jU^)=jRiz)~I. (11.1) 

Then by mM . 

j^{z) = 0{N~'^) uniformly on the circles, 

j'^{z) = 0{e^'^^^^^) for some c{z) > 0, on the remaining parts of Tr, 
where c(x) satisfies (16. 5p as x — oo. We can apply the following general result. 
Proposition 11.1. Assume that v{z), z E Tr, solves the equation 

= ^(!MM,„, ,erH, (11.3) 

2m r„ z_ — u 



R 



;ii.2) 



where z_ means the value of the integral on the minus side o/Tr. Then 

R[^)^i-1-! hMzM,„, ,ec\r«, (11.4) 

2711 Jy^ z-u 

solves the following RH problem: 

(i) R{z) is analytic on C\Tr. 

(ii) R+{z) = R.{z)jR{z),zeTR. 

(iii) R{z) = I + 0{z-^) , z ^ oo. 

Proof. From ffTL3|) . ffTTD . 

R4z) = v{z), zeTR. (11.5) 
By the jump property of the Cauchy transform, 

R^z) - R.{z) = v{z)jI{z) = R-{z)jI{z), (11.6) 

hence R+{z) = R-{z)jR{z). From (111.40 . R{z) = I + 0{z^^). Proposition II 1. II is proved. 
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Equation (111.31) can be solved by perturbation theory, so that 

oo 

t;(^) = J + (11.7) 

A;=l 

where for > 1, 



- ) 



/ !!tl(!MM<i„. ,gr«. (11.8) 



R 



and Vq{z) = I. Series (I11.7P is estimated from above by a convergent geometric series, so it 
is absolutely convergent. From (lll.2p we obtain that there exists C > such that 

Observe that 

vi{z) = -^[ ^^du, zeTji. (11.10) 
We apply this solution to find Rn{z). The function Rn{z) is given then as 

oo 

Rn{z) =1 + J2RNk{z), (11.11) 



where 



In particular. 



/ "J^^MlMdu. (11.12) 



R 



Rm{z) = -^. [ ^-^du. (11.13) 
2m Jy^^z-u 



From flll.yp we obtain that there exists Co > such that 



Hence from (111.111) we obtain that there exists Ci > such that for A; > 0, 

N''+\l + \z\) 



Rn{z) = I + ^RNj{z) +eNk{z), \eNk{z)\ < ^^^-^^^ ^ i^n " (11.15) 



In particular. 



R^{z)=I + 0[j;^^^-^] asiV^oo, (11.16) 



uniformly for 2; G C \ Fij. 
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12. ASYMPTOTICS OF THE RECURRENT COEFFICIENTS 



Let us summarize the large N asymptotics of orthogonal polynomials. From (110. ip and 
(111.161) we obtain that 



Sn{z)= [I + O 
Sn(z)= {I + O 



N{\z\ + 1) 



M{z), zeC\D, 
Un{z), zeD; 



;i2.i) 



D = [j[Dia„r)UDib„r)]. 



From (17. 6p we have that 



Sn{z) y^-NG{z) ij, ze = U>C+, 



Tm(z) 



Sn(z) 



1 



ec- = [jq 



;i2.2) 



Finally, from (16.21) we obtain that 

1 



Ym(z) = ( 



e 2 



e 2 



■0-3 



■0"3 



e 2 



I + O 
I + O 
I + O 



Ni\z\ + 1] 



N{\z\ + 1) 



N{\z\ + 1) 



M{z) 



1 



\ zeC^\D, 



M{z)e''\'^'^-^y\ z eC\{DUC+UC~). 



(12.3) 

This gives the large N asymptotics of the orthogonal polynomials and their adjoint functions 
on the complex plane. Formulae ( 13.17^ and ( 13.18^ give then the large N asymptotics of the 
recurrent coefficients. Let us consider 7^. 
From (16.21) we obtain that for large z, 



/ + — + 4 + 

z z^ 



hence 
and 



= Yn{z)z-''''^ = e^'''T{z)e^i''^'^-^-'"^'h 

= ef -3 (j^Il^Tl^ \ ^N[giz)-l-losz]a,^ 

\ z z"^ J ' 



ll=[Yl]l2[Yl]21 = [Tl]l2[Tl]21. 

From ffT^ . ffT^ we obtain further that 

7^ = [Mi]i2[Mi]2i + 0(iV-i), 



;i2.4) 

;i2.5) 
;i2.6) 

;i2.7) 
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and from f l8.5ip . 



[Mi]i2 



6{—u{oo) + d + c)6{u{oo) + d) {bj — aj) 
9{u{oo) + d + c)9{-u{oo) + d) j-^ (-4i) 

6{u{oo) — d + c)6{—u{oo) — d) {bj — aj) 
9{-u{oo) - d + c)9{u{oo) - d) Ai 



(12.8) 



hence 



7^ 



4 



EC-. 



■1 «jJ 



^2(^(00) + d)d{u{oo) + ^ - d)d{-u{oo) + ^ + d) 

e\u{oo) - d)e{-u{oo) + f^- d)e{u{oo) + f^ + d) 



+ 0{N~'] 



where d is defined in (18.451) . Consider now Pn-i- 
From (112. 4p we obtain that 



[Villi = [Tilii + Ng,, [Y^hi = e-''\[T2]2i + [TihiNg^ 



hence 



N-l 



[^2)21 r^n _ [^2] 21 ^rp 1 

- i^ljll — 77771 [-tljll) 



[Yi] 



21 



and by ([m]), (^T^, 
From flHIlH]) we find that 



Af-l 



\M2\2i 
[Mi]2l 



[T^l]21 

~[Mi]n + 0{N~'). 



(u(cxo) — d) 



+ 



9{u-d + c) 



e{u - d) 



,u (00) 



u=u{co) 



[Ml] 



11 



9{u{oo) + d) 
9{u{oo) + d + c) 



\/e(u(oo) - d + c) Ve(u(oo) - d + c) , 
9{u{oo)-d + c) 9{u{oo) - d + c) ' ^ ' 

6{u + d + c) 



e{u + d) 



,u [00] 



u=u{oo) 



Hence, 



P. 



N-l 



Ve{u{oo) + d + c) Ve(u(oo) + d + c) . 
9{u{oo) + d + c) 9{u{oo) + d + c) ' ^ ' 



E,'=i(&?-S') , / Ve(Moo) + - rf) V^^(n(oo) + + rf) 



u{oo) + ^-d) e{u{oc) + ^ + d) 



6{u{oo) + d) 6{u{oo) - d) 



(12.9) 



(12.10) 



(12.11) 



(12.12) 



(12.13) 



(12.14) 
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This formula can be also written in the shorter form, 



/9, 



N-l 



+ 



log- 



u{z) + - d)e{u{z) + d) 



e{u{z) + + d)e{u{z) - d) 



0{N 



In the one-cut case, g = 1, ai = a, 61 = 6, formulae fll2.9p . (112. 14p simplify to 

a + h 



+ 0{N~^). 



(12.15) 



;i2.16) 



Formula (112. 9p is obtained in [DKMVZ2]. Formula (112. 14p slightly differs from the formula 
for in [DKMVZ2]: the first term, including a^-'s, 6j's, is missing in [DKMVZ2]. 

13. Universality in the random matrix model 



By applying asymptotics (112. 3p to reproducing kernel (I3.25P , we obtain the asymptotics of 
the eigenvalue correlation functions. First we consider the eigenvalue correlation functions 
in the bulk of the spectrum. Let us fix a point Xq e Int J = Uj^]^(aj, bj). Then the density 
Pvi^o) > 0. We have the following universal scaling limit of the reproducing kernel at Xq: 

Theorem 13.1. .4s ^ 00, 

1- ^ f ^ \ sin[7r(M — i;)] /,o.n 

hm — —Kn{xo + — — ^,Xo + ^^--^-^] = — (13.1) 



N^oc Npvixo) V Npvixo) Npvixo)J n{u - v) 

Proof. Assume that for some I < j < q and for some e > 0, we have {xq, x, y} G {aj+e, bj—e). 
By (13:251) and dOj) . 



_ NV(x) _ NV(y) 

Kn{x, y) = \_],^ (0 1) Yj,liy)Y^^ix) 



2TTi{x — y) 

_ NV(x) _ NVjy) 
6 2 6 2 



e^K(^)-^])r^;(z/)T^+(x) 



,N[9+(x)- 





(13.2) 



2m[x — y) 

Now, from (17.60 we obtain that 

NV(x) NV(y) 

KNix,y) = " (0 e^K(^)"^]) My)S],liy)S^+ix)j-\x) ^ ^ 



^N[g+{x)-L] 



NV(x) NV{y) 
6 2 6 2 



2m{x — y) 



^N[-G{y)+g+{y)-L] ^N[g+{y)-L]\ 5^^^) 



X Sn+{x) 



^N[g+(x)-L] 
,iV[G(x)+g+(x)-i; 



By dEISD, 



V{x) . , / 
4^ + ^-(^)-2^^ 



G{x) 



hence 



KN{x,y) 



1 



27ri{x — y) 



NG(y) 

6 2 



NG(y) 



S^\{y)SN+{x) 



NG(x) 

6 2 

NGjx) 
— 6 2 



(13.3) 
(13.4) 

(13.5) 
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By (imlD, 

Sjv+ix) = Rn{x)M+{x). (13.6) 

Observe that M+(a;) and Rn{x) are analytic on {aj + e, bj — e) and Rn{x) satisfies estimate 
(111.161) . This imphes that a.s x — y ^ 0, 

S],l{y)SMx)=I + 0{x-y), (13.7) 

uniformly in A^. Since the function G{x) is pure imaginary for real x, we obtain from (113.51) 
and (nJl\\ that 

K.ix, y) = [e^^ - e^^J +0(1). (13.8) 

By dEH, 

N[G{x) - G{y)] 



mN pv{s)ds = 7rtNpv{0{x-y), ^ e [x,y], (13.9) 
^ Jy 

hence 

sin[7rArpv'(0(a^-2/)] .^/.^ nqin^ 
/s:^(x,2/) = + 0(1). (13.10) 

7c{x - y) 



Let 



u V 
^ = ^o + ^j — y = XQ + — — -— , 13.11 



where u and v are bounded. Then 



1 "'I +0(iV-). (13,12) 



A''py(xo) ' 7r(M - f ) 

which implies (113.11) . □ 

Consider now the scaling limit at an edge point. Since the density py is zero at the 
edge point, we have to expect a different scaling of the eigenvalues. We have the following 
universal scaling limit of the reproducing kernel at the edge point: 

Theorem 13.2. If Xq = bj for some 1 < j < q, then for some c > 0, as N ^ oo, 

^ r. f u V \ Ai{u)Ai'{v) - Ai'{u)Ai{v) 

Similarly, if xq = aj for some I < j < q, then for some c> 0, as N oo, 

I u V \ Ai(u)Ai'(v) - Ai'(u)Ai(v) 

The proof is similar to the proof of Theorem 113.11 and we leave it to the reader. 

Lecture 3. Double scaling limit in a random matrix model 
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14. Ansatz of the double scaling limit 

This lecture is based on the paper [BI2]. We consider the double- well quartic matrix 
model, 

Aiiv(M) = Z-ie-^^^(*^)dM (14.1) 

(unitary ensemble), with 

V{M) = ^ + —, t<0. (14.2) 

The critical point is = —2, and the equilibrium measure is one-cut for t > —2 and two-cut 
for t <-2, see FigJl 

The corresponding monic orthogonal polynomials Pn{z) = + . . . satisfy the orthogo- 
nality condition, 

/oo 

P„(z)P„(z)e-^^(^)dz = hj^^. (14.3) 
-oo 

and the recurrence relation, 

zPniz) = Pn+l{z) + RnPn^l{z). (14.4) 

The string equation has the form, 

Tl 

Rnit + i?„_i + i?„ + Rn+l) = — . (14.5) 

For any fixed e > 0, the recurrent coefficients Rn have the scaling asymptotics as ^ oo: 

Rn = aQ)+{-irbQ)+0{N-'), e<^<X^-e, (14.6) 

and 

Rn = a[^)+OiN-'), s-'>^>X, + e, (14.7) 

where 

4 • 

The scaling functions are: 



(14.8) 



a(A) = -| bW = '^^\ A<Ae, (14.9) 



and 



a(A) = + V^^ . A>A.. (14.10) 

Our goal is to obtain the large asymptotics of the recurrent coefficients Rn, when n/N 
is near the critical value Ac. At this point we will assume that t is an arbitrary (bounded) 
negative number. In the end we will be interested in the case when t close to (—2). Let us 
give first some heuristic arguments for the critical asymptotics of Rn- 
We consider A — > oo with the following scaling behavior of ^: 

^ = K + c,N~^l\ Co=(^)'^', (14.11) 
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where y G (— C)0, oo) is a parameter. This hmit is called the double scaling limit. We make 
the following Ansatz of the double scaling limit of the recurrent coefficient: 

Rn = -\ + N~^l\_iY^^^^y^ ^ N-''/\^v{y) + O(iV-i), (14.12) 

where 

c, = mf\ ^2 = ^^) . (14.13) 

The choice of the constants co,ci,C2 secures that when we substitute this Ansatz into the 
left-hand side of string equation (114.51) . we obtain that 

/?„(t + /2„„i + i?„ + i?„+i) = ^ + Ar-2/3co [v - 2u' - y) +N-\-lY {u" -uv) + ... (14.14) 

By equating the coefficients at N^'^/^ and A^^^ to 0, we arrive at the equations, 

v = y + 2u^ (14.15) 

and 

u" = yu + 2u^, (14.16) 

the Painleve II equation. The gluing of double scaling asymptotics (114. 12p with (114. 6p and 
(I14.7P suggests the boundary conditions: 

u ~ C\J—y^ y ^ —oo; u ^ 0, y oo. (14.17) 

This selects uniquely the critical, Hastings-McLeod solution to the Painleve II equation. 
Thus, in Ansatz (114.121) u{y) is the Hastings-McLeod solution to Painleve II and v{y) is 
given by (114.151) . The central question is how to prove Ansatz (114. 12p . This will be done 
with the help of the Riemann-Hilbert approach. 

We consider the functions ipniz), n = 0,1, . . ., defined in (12.261) . and their adjoint functions, 

2tc u- z 

We define the Psi-matrix as 

v^„,(.)=f/^"(f. ^"(f.) (14.19) 

The Psi-matrix solves the Lax pair equations: 

K{z) = NA^{z)^^{z), (14.20) 

^r.+l{z) =Un{z)^n{z). (14.21) 

In the case under consideration, the matrix is given by formula (12.580 . with g = 1: 

An{z) = (~ (f + * + "^^O ^"^'^f ^^A, e^ = t + R^ + R^^,. (14.22) 
\-R'J\z^ + en-i) 'i + ^ + zRj 

Observe that (114. 20p is a system of two differential equations of the first order. It can be 
reduced to the Schrodinger equation, 

- 7]" + N^Ut] = 0, 7] = , (14.23) 

vail 
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where the matrix elements of An{z), and 



U = -detAn + N-^ 



an. 



air 



ai2 



N 



-2 



jau)" 3((ai2) 
2ai2 4(012) 



/\2-\ 



14.24) 



see [BIl], [BI2]. 

The function \l/„(-2) solves the following RHP: 

(1) \l'n(-2) is analytic on {Imz > 0} and on {Im^; < 0} (two-valued on {Im^; = 0}). 

, X G 



(2) ^„+(a;) = ^„_(a:) 

(3) As z — > 00, 







2; ^ 00, 



where Fz,, k = 0,1, 2, 



are some constant 2x2 matrices, with 







R, 



-1/2 



Ti 



1 

R'J' 



(14.25) 



14.26) 



A. 



In h„ 



2 ' 



and 



1 
-1 



is the Pauli matrix, 



Observe that the RHP implies that det\l/„(z) is an entire function such that det\l/„(oo 

det^„(z) = i?;^/^ zeC. 



Rri^^'^ , hence 



(14.27) 

We will construct a parametrix, an approximate solution to the RHP. To that end we use 
equation (114.201) . We substitute Ansatz (114.121) into the matrix elements of A^ and we solve 
(114.201) in the semi classical approximation, as — 00. First we determine the turning 
points, the zeros of detAn{z). From (114.221) we obtain that 



deiAn{z) = an{z) = - 
Ansatz flM.lip . 014.121) implies that 
^ - A. = coN-^/'y 

hence 

det A„ (z) - 



z 

T 



(^J^ — Acj z"^ + RnOnOn+l, 



tz'' 



en = 2c2N-'/My) + o{N-'). 



(14.28) 



(14.29) 



(14.30) 



^ ^ + coN-'/^yz' - 2t[c2v{y)]^N-^/^ + 0(Ar-5/3). 

We see from this formula that there are 4 zeros of det An, approaching the origin as A^ — > 00, 
and 2 zeros, approaching the points ±zo, zq = ^ —It. Accordingly, we partition the complex 
plane into 4 domains: 

(1) a neighborhood of the origin, the critical domain , 

(2) 2 neighborhoods of the simple turning points, the turning point domains fij^f , 

(3) the rest of the complex plane, the WKB domain 

We furthermore partition into three domains: ^^If^ and see Fig|8l 
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,WKB 
CO 




Figure 8. The partition of the complex plane. 

15. Construction of the parametrix in 
In fi^^^ we define the parametrix by the formula, 



(15.1) 



where Cq ^ 0, Ci are some constants (parameters of the solution). To introduce T{z) and 
^{z), we need some notations. We set 



as an approximation to Rn, and 



AVz) 



6^=t + RI + 



as an approximation to An{z). We set 



<(.) = - ^ + - A.) + N~'/^-ty/'2-'/My? - My? 

- N-'/%-in-2t)-'/My)], Hv) = n{y), 
as an approximation to det A„(z). Finally, we set 

[/0 = -a°(;2) + iV-i 



('^12)^ 

11 
^12 



as an approximation to the potential U in (114.241) . With these notations, 

= r fi{u) du, ii{z) = ^/uKz). 



(15.2) 



(15.3) 



(15.4) 



(15.5) 



(15.6) 



2]V 
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where is the zero of U^{z) which approaches zq as N oo. Also 

n^) = {^Y'( h detT(z) = l. (15.7) 



From (115. ip we obtain the following asymptotics as 2; ^ oo: 
where 



A? = lim 



Na.)-{^-nlnz 



(15.8) 
(15.9) 



The existence of the latter limit follows from (115. 6p . 
In the domains fl^^^ we define 

^WKB^^) ^ ^WKB^^)^^^ ±Im2>0, (15.10) 

where \l/^^^(2;) is the analytic continuation of \l/^^^(z) from ^^^^^ to Q^^^, from the 
upper half-plane, and 

5.^(1 7), s_^(;;). ,15.11) 

Observe that \l/^™(2;) has jumps: 

V^WKB^^) ^^j^ 0(e-=^))vl/WKB(^)^ ^ g Q^WKB ^ ^g^WKB ^ g^WKB^^ ^ ^^^^^^^ 

and 

Vl/WKB^^) ^ Vl>WKB(^) (0 T) ' ^ ^ ^ ^ ^ • (^^-^2) 

16. Construction of the parametrix near the turning points 
In fij^ we define the parametrix with the help of the Airy matrix-valued functions, 

ii.(--) = (!?S!! S'^i!!), (16.1) 



where 

Vq^z) = Ai{z), yi{z) = e-'^/^Ai {e-^'^'/^z) , y2{z) = e'^/^Ai {e^'^'/^z) . (16.2) 

Let us remind that Ai(2;) is a solution to the Airy equation y" = zy, which has the following 
asymptotics as 2; — > 00: 

1 / 2z^/^ \ 
Ai(^) = exp + 0(|2|-^/2) , -TT + e <aTgz <7r -e. (16.3) 



2v/?zi/4 y 3 
The functions yj{z) satisfy the relation 

yiiz) - y2{z) = -iyoiz). (16.4) 
We define the parametrix in fij^ by the formula, 

^/^^{z) = W{z)N'"^^Yi^2{N'^^^w{z)), ±lmz>0, (16.5) 
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where 



2/3 



5 



(16.6) 



with C,{z) defined in fll5.6l) above. Observe that w{z) is analytic in The matrix- valued 
function W{z) is also analytic in fij^, and it is found from the following condition of the 
matching ^^^{z) to ^wkbJ^) qqtp. 

see [BIl], [BI2]. A similar construction of the parametrix is used in the domain Qj^ . 

17. Construction of the parametrix near the critical point 

17.1. Model solution. The crucial question is, what should be an Ansatz for the parametrix 
in the critical domain To answer this question, let us construct a normal form of system 

of differential equations (114. 20p at the origin. If we substitute Ansatz (114.121) into the matrix 
elements of An{z), change 

^{z) = <l>{CN^/h), C=^-^, (17.1) 

and keep the leading terms, as ^ oo, then we obtain the model equation (normal form), 

<l>'(s) = A(s)<l>(s), (17.2) 

where 

/ {-irAu{y)s 4s' + {-ir2w{y) + v{y)\ 

^^''^ ~ \-As^ + {-ir2w{y)-v{y)) -{~l)"Au{y)s J' ^^''^^ 

and w{y) = u'{y). In fact, this is one of the equations of the Lax pair for the Hastings-Mcleod 
solution to Painleve II. Equation (117. 2p possesses three special solutions, $j, j = 0, 1, 2, which 
are characterized by their asymptotics as |s| — > oo: 



f^Hs)\ f cos i^s^ + ys \ , , ^ vr 



\n+l 



(17.4) 



The functions $^'^(s) are real for real s and 

$i(-s) = (-l)"$i(s), $2(-s) = -(-l)"<l>2(s). (17.5) 
We define the 2x2 matrix-valued function on C, 

$(5) = ($o(s), $1,2(5)) , ±Ims > 0, (17.6) 

which is two-valued on M, and 



$+(s) = $_(.)(i /), seR. (17.7) 



^0 1 

The Ansatz for the parametrix in the critical domain is 

vI/CP(^) = CoViz)^ {N'/'Ciz)) , (17.8) 
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where Cq is a constant, a parameter of the solution, ({z) is an analytic scalar function such 
that C'(-2) 7^ and V{z) is an analytic matrix-valued function. We now describe a procedure 
of choosing ({z) and V{z). The essence of the RH approach is that we don't need to justify 
this procedure. We need only that ({z) and V{z) are analytic in and that on dQ'"^ , 
Ansatz ffTTISD fits to the WKB Ansatz. 

17.2. Construction of C(z): Step 1. To find ((z) and V(z). let us substitute Ansatz (117.81) 
into equation f l2.54p . This gives 

V{z) [C{z)N-^/^A {N'/\{z))] V-\z) = An{z) - N-'V'{z)V~\z) (17.9) 

Let us drop the term of the order of A^^^ on the right: 

V{z) [az)N-'/'A {N'/'dz))] V'\z) = A^{z), (17.10) 

and take the determinant of the both sides. This gives an equation on ( only, 

[C'W]V(CW) = a„(^), (17.11) 

where 

/(() = det A{N^/^C) = 16C' + 8Ar-2/3y(2 ^ N~i/3^y2^y^ _ 4^2^^^] ^^^^^2) 

and 

tz^ Z^ / Tl \ 

ttniz) = det ArXz) = -— -^+[j^->'c)z^ + (17.13) 

where 



en=t + Rn + Rn+1. (17.14) 

Equation (114. 12p implies that 

On = 2c2N-^'^v{y) + O(Ar-i). 

At this stage we drop all terms of the order of A^~^, and, therefore, we can simplify / and 
On to 

/(C) = 16C' + 8Ar-2/V' (17.16) 

and 

an{z) = --^-'^+{^-^) (17.17) 

Equation (117.111) is separable and we are looking for an analytic solution. To construct an 
analytic solution, let us make the change of variables, 

z = CN-^'^s, C = N-^'^a. (17.18) 

Then equation (I17.1ip becomes 

[a'{s)fU<y{s)) = ao(.), (17.19) 

where 

Ma) = 16a^ + 8ya^, ao{s) = 16s^ + Sc^^N^/^ - A,) - N-^'^cs\ (17.20) 

When we substitute equation (114. 13p for y, we obtain that 

ao(s) = 16/ + Sys^ - N~'^/^cs^. (17.21) 
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When y = 0, equation (117. 19p is easy to solve: by taking the square root of the both sides, 
we obtain that 

1/2 



2 / 

a a = s 



2(1- liV-2/3cs2 j , (17.22) 



hence 



1 1/3 

(17.23) 



is an analytic solution to (117.191) in the disk |s| < eN^^^, for some e > 0. This gives an 
analytic solution C(^) = N-'^/^a{C-^N^^^z) to equation fll7.13|) in the disk \z\ < Ce. 

When y ^ 0, the situation is more complicated, and in fact, equation (117.191) has no 
analytic solution in the disk |s| < eN^^^. Consider, for instance, y > 0. By taking the 
square root of the both sides of (117.191) . we obtain that 



a a 



.2 , yV^', 



2/ V 2 16 



cr' = s(s' + ^- i-N~^l^cs^ , (17.24) 



1/2 



The left hand side has simple zeros at icrg = iiy^, and the right hand side has simple 
zeros at ±so, where Sq = o"o + 0{N^'^^^). The necessary and sufficient condition for the 
existence of an analytic solution to equation (I17.19P in the disk \s\ < eN^^^ is the equality 
of the periods. 



^0 / „ 1 \ 1/2 



Pi= I ^ f + da = P2= r + ^N-^'^cs^ ] ds, (17.25) 



.0-27 V 2 16 

and, in fact. Pi ^ P2- To make the periods equal, we slightly change equation (114.131) as 
follows: 

y = Co ^iV^/^ - A,) + a, (17.26) 

where a is a parameter. Then 

/•■50(a) / „_ 1 \ 

P2 = PM) ^ sls' + ^- T^N~"^'^n ds. (17.27) 

J-so{a) V 2 16 J 

It is easy to check that P2{0) 7^ 0, and therefore, there exists an a = 0(A^-2/3) such that 
Pi = P2. This gives an analytic solution cr(s), and hence an analytic C{z). 

17.3. Construction of V{z). Next, we find a matrix- valued function V{z) from equation 
(117.101) . Both V{z) and V~^{z) should be analytic at the origin. We have the following 
lemma. 

Lemma 17.1. Let B = {bij) and D = {dij) he two 2x2 matrices such that 

Tr5 = TrD = 0, det5 = detL'. (17.28) 
Then the equation VB = DV has the following two explicit solutions: 
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We would like to apply Lemma [17.11 to 

B = C'(z)iV-2/3^ {N^/^C{z)) , D = Ar,{z). (17.30) 

The problem is that we need an analytic matrix valued function V{z) which is invertible in 
a fixed neighborhood of the origin, but neither Vi nor V2 are invertible there. Nevertheless, 
we can find a linear combination of Vi and V2 (plus some negligibly small terms) which is 
analytic and invertible. Namely, we take 

Viz) = rr^r^ Wiz) (17.31) 
A/det W{zj 

where 

Wiz) = (J^'f'^ ->/^) - ^}'^'\- Yf\~ , (17.32) 

\bu{z) ~ dii{z) - a2iz, 612(2;) - ^21(2) - ^22 / 

and the numbers aij = 0{N^^) are chosen in such a way that the matrix elements of W 
vanish at the same points ±Zq, zq = 0{N~^^^), on the complex plane. Then V{z) is analytic 
in a disk \z\ < e, e > 0. 

17.4. Construction of C{z): Step 2. The accuracy of C(^)) which is obtained from equation 
(117.111) . is not sufficient for the required fit on |z| = e, of Ansatz (117.81) to the WKB Ansatz. 
Therefore, we correct ({z) slightly by taking into account the term —N~^V'{z)V~^{z) in 
equation (117.91) . We have to solve the equation, 

[C{z)]^N~^/^ det A (iV^/'C(^)) = det [An{z) - N-W{z)V~\z)] . (17.33) 

By change of variables (117. 18p , it reduces to 

[a'{s)rh{a{s)) = a,{s), (17.34) 

where 

f,(a) = I6a' + 8ya' + [v\y) - Aw\y)]; 

a,{s) = 16s^ + 8iy-a)s^ + [v\y)-4w\y)]+rM{s), rM{s) = 0{N~^/^). 

The function /i(cr) has 4 zeros, icXj j = 1,2. The function ai(s) is a small perturbation of 
/i(s), and it has 4 zeros, ±Sj, such that \sj — aj\ — >• as iV — >• cxd. Equation (I17.34p has an 
analytic solution in the disk of radius eN'^^^ if and only if the periods are equal, 

/(Jf) t-Sj 
y^W)da = P2,= / ^^)ds, J = 1,2. (17.36) 
-CTj J-Sj 

To secure the equality of periods we include ai{s) into the 2-parameter family of functions, 
ai(s) = 16s^ + 8(1/ - a)s^ + [v^{y) - 4:W^{y)] + tn^s) + /?, (17.37) 
where —00 < a,/? < 00 are the parameters. A direct calculation gives that 

[ dP21 dP2l \ 

det ^ U 0, (17.38) 

\ da dp J 

see [BI2], hence, by the implicit function theorem, there exist a,/? = 0(A^~^/^), which solve 
equations of periods (117.361) . This gives an analytic function C{.z), and hence, by (117.81) . an 
analytic '^^^{z). 
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Figure 9. The contour T. 

The function "^^^(z) matches the WKB-parametrix "^^^^{z) on the boundary of fl^^ . 
Namely, we have the following lemma. 

Lemma 17.2. (See [BI2].) // we take Cq = Cq and Ci = -| lni?° then 

^CP(^) = (J + 0(Ar^i))^™(z), z G Sr^^P. (17.39) 

We omit the proof of the lemma, because it is rather straightforward, although technical. 
We refer the reader to the paper [BI2] for the details. 

17.5. Proof of the double scaling asymptotics. Let us summarize the construction of 
the parametrix in different domains. We define the parametrix \E'° as 



WKB y qWKB 



(17.40) 



vi/TP(^), zen^unr, 

where ^wkbj^^^ gj^^^ f[T51|l . fITCTD . <i/^^{z) by ffT63D . and by fHTS]) . Consider 
the quotient, 

X(^) = vI/„,(^)[xl/0(^)]-i. 
^(z) has jumps on the contour F, depicted on Figj9l such that 



X+{z) = X_{z) [I + oiN-\i + \z\y 

From (114.251) and (115.81) we obtain that 



(17.41) 
(17.42) 



Xiz) =Xo + 



0(z~ 



oo, 



where 



V2 



and 



1 



p^g{Ci+Aa-A„)a3 _ roe('^i+^""^")'^='(i?°)^/^ 



0-1 



0-1 



1 

1 



;i7.43) 
17.44) 

;i7.45) 



The RHP shares a remarkable property of well-posedness, see, e.g. [BDT], [CG], [LiS], [Zh]. 
Namely, equations 017.421) . fll7.43p imply that 



X^^X{z)=I + 0{N~\l + \z\)-^), ze 



;i7.46) 
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This in turn implies that 
or, equivalent ly, 



(c,+A0-A„).3r-ir^e(c,+A0-A„)<x3 _ (i?0)i/2^^ ^ o(iv-i). 



By (fTCTIl . 



i Q i le 



Q g-2(Ci+A0-A„) 
p ^2(Ci+A0-A„) 



hence fll7.48p reduces to the system of equations, 

^^g2(C,+A0-A„)^ (^0)1/2 

By multiplying these equations, we obtain that 

/2„ = i?°+0(iV-i). 
This proves Ansatz (114. 12p . Since Ci = — | lni?°, we obtain from (117.511) that 



.2(A„-A0) 



1 + OiN-^) 



or equivalently, 



where 



hr, = exp (^2N ^i{u) du^ (l + 0{N-^)) , 

V{z) nlnz 



jj,{u) du = lim 



2]V 



N 



Thus, we have the following theorem. 



(17.47) 
(17.48) 

(17.49) 

(17.50) 

(17.51) 
(17.52) 
(17.53) 
(17.54) 



Theorem 17.3. (See [BI2]). The recurrent coefficient Rn under the scaling ( 14- 11 ) has 
asymptotics \14-li^ - The normalizing coefficient has asymptotics ([777 



zeC. 



(17.55) 



Equations 017.411) and 017.461) imply that 

^„(z) =Xo[l + 0{N-\1 + \z\)-'] 

The number Co is a free parameter. Let us take Cq = 1. From 017.44p and 017.5ip we obtain 
that 

(^0)1/4 



V2 



{1 + 0{N-')) 



hence 



I' DO \ 1/4 



zee. 



;i7.56) 



;i7.57) 



This gives the large N asymptotics of ^^.^(-z) under scaling 014.110 . as well as the asymptotics 
of the correlation functions. In particular, the asymptotics near the origin is described as 
follows. 
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Theorem 17.4. (See [BI2]). Let $o(;z;y) = {^2(^^^.y^j solution for n = to system 



(17.2), with the asymptotics at infinity as in ( 17.4\ )- Then the following double scaling limit 



holds: 



det {Qc{ui,Uj]y))^.-^ 



,j=l,...,m 



where c = C'(0) > 0, and 

ftK = (17.59) 

7r{u — V) 

Let us mention here some further developments of Theorems 1 1 7. 31 [1 7. 4[ They are extended 
to a general interaction V{M) in the paper [CK] of Claeys and Kuijlaars. The double 
scaling limit of the random matrix ensemble of the form Z]^^\ det M\'^°'e~^'^^^^'^^^dM, where 
a > — I , is considered in the papers [CKV] of Claeys, Kuijlaars, and Vahlessen, and [IKO] 
of Its, Kuijlaars, and Ostensson. In this case the double scaling limit is described in terms 
of a critical solution to the general Painleve II equation q" = sq + 2q^ — a. The papers, 
[CK], [CKV], and [IKO] use the RH approach and the Deift-Zhou steepest descent method, 
discussed in Lecture 2 above. The double scaling limit of higher order and Painleve II 
hierachies is studied in the papers [PeS], [BE], and others. There are many physical papers 
on the double scaling limit related to the Painleve I equation, see e.g., [BrK], [DS], [CM], 
[DGZ], [Wit], and others. A rigorous RH approach to the Painleve I double scaling limit 
is initiated in the paper [FIK] of Fokas, Its, and Kitaev. It is continued in the recent 
paper of Duits and Kuijlaars [DuK], who develop the RH approach and the Deift-Zhou 
steepest descent method to orthogonal polynomials on contours in complex plane with the 

exponential quartic weight, e~^^^^^\ where t < 0. Their results cover both the one-cut 
case —j^<t<0 and the Painleve I double scaling limit at t = — . 

Lecture 4. Large N asymptotics of the partition function of random 
matrix models 

18. Partition function 
The central object of our analysis is the partition function of a random matrix model, 

/oo I'oo N—l 

■■■ n (-^^^ - ^fc)'e-'^^^"-^(^^^rf^i ...dzN = N\l[ hn, (18.1) 

■oo J -co i<j^k<N n=0 

where V{z) is a polynomial, 

2d 

V{z) = Y,^,z^, V2d>0, (18.2) 
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and hn are the normalization constants of the orthogonal polynomials on the line with respect 
to the weight e~^^^^\ 

/oo 

P„(z)P„(^)e-^^(^)rfz = hJnm, (18.3) 
-oo 

where Pn{z) = + . . . We will be interested in the asymptotic expansion of the free energy, 

F^ = -^\nZ^, (18.4) 

as ^ oo. 

Our approach will be based on the deformation of V{z) to z"^, 

Tt : V{z) V{z; t) = {1 - t-^)z^ + V{t-^'\), (18.5) 

1 < t < oo, so that 

nV{z) = V{z), T^V{z) = z\ (18.6) 

Observe that 

TtTs = Tts- (18.7) 

We start with the following proposition, which describe the deformation equations for 
and the recurrent coefficients of orthogonal polynomials under the change of the coefficients 
ofV{z). 

Proposition 18.1. We have that 

]_d\nhn_ _ 

N dvk 

^ = ? ([Qln-l,n-l - mnn) , (18.8) 

/n[^; \n,n~l /n+1 \n+l,ni 

where Q is the Jacobi matrix defined in Ii2.33\) . 

The proof of Proposition 1 1 8 . 1 1 is given in [BI3]. It uses some results of the works of Eynard 
and Bertola, Eynard, Harnad [BEH]. For even V, it was proved earlier by Fokas, Its, Kitaev 
[FIK]. 

We will be especially interested in the derivatives with respect to V2- For k = 2, Proposition 
118.11 gives that 

1 dlnhn 2 n2 2 



N dv2 



In Pn In+li 



^|^ = y(7Li+/5Li-7^.i-/?^), (18.9) 

1 djSn 2 n , 2 n 2 n 2 n 

iV 'dv^ " 7„P„-1 + 7nPn - In+lPn " 7n+lPn+l- 

Next we describe the f2-deformation of Z^- 

Proposition 18.2. (See [BI3].) We have the following relation: 

Jf2 = 7^(7^1 + 7^+1 + Pl + WnPn-1 + /?^-i). (18.10) 
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Observe that the equation is local in A^. For V{z) = V2Z^ + ^42;^, Proposition 118.21 was 
obtained earlier by Its, Fokas, Kitaev [FIKO]. Proposition 118. 21 can be applied to deformation 
fll8.5p . Let 7n(T), (inij^ be the recurrence coefficients for orthogonal polynomials with respect 
to the weight e~^^(^''^\ and let 

/OO /"OO 
•••/ n (%--Sfc)V^^^"^^'^'^^-2l---^^-2jV (18.11) 

■00 J-cxD x<j<k<N 

be the partition function for the Gaussian ensemble. It can be evaluated explicitly, and the 
corresponding free energy has the form, 

.Gauss ^_^.J (2vr)^/^ 
iV2 



^'''((i^n«y (18.12) 

By integrating twice formula (118. lOp . we obtain the following result: 
Proposition 18.3. 

= F^— + / {7^(r) [7^_i(r) + 7^+1 (r) + /3^(r) + 2/?;v(r)/3^_i(r) 

(18.13) 

19. Analyticity of the free energy for regular V 

The basic question of statistical physics is the existence of the free energy in the thermody- 
namic limit and the analyticity of the limiting free energy with respect to various parameters. 
The values of the parameters at which the free energy is not analytic are the critical points. 
When applied to the "gas" of eigenvalues, this question refers to the existence of the limit, 

F = lim Fm, (19.1) 

and to the analyticity of F with respect to the coefficients Vj of the polynomial V. The 
existence of limit (119.10 is proven under general conditions on V, not only polynomials, 
see the work of Johansson [Joh] and references therein. In fact, F is the energy Ey of 
minimization problem (14. 7p . (14. 8p . so that 

F = Iviiyv), (19.2) 

where uy is the equilibrium measure. The following theorem establishes the analyticity of 
F for regular V. We call V regular, if the corresponding equilibrium measure uy is regular, 
as defined in (I4.17p . (14.181) . We call V, g-cut regular, it the measure is regular and its 
support consists of q intervals. 

Theorem 19.1. (See [BI3].) Suppose that V{z) is a q-cut regular polynomial of degree 2d. 
Then for any p < 2d, there exists ti > such that for any t G [— ti, ti], 

(1) the polynomial V{z) + tz^ is q-cut regular. 

(2) The end-points of the support intervals, [aj(t), 6j(t)], z = l,...,g, of the equilibrium 
measure for \y{z) + tz"^] are analytic in t. 

(3) The free energy F{t) is analytic in t. 
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Proof. Consider for j = 0, . . . , q, the quantities 

T,(a, b;t) = ^(f ^'^'r^^ dz, R{z) = ([{z - a,){z - 6.), (19.3) 



z] 



i=l 



where F is a contour around [ai, h^. Consider also for = 1, . . . , g — 1, the quantities 

Nk{a,h-t) = 7T- { h{z-t)^R{z)dz, (19.4) 

where Vk is a contour around [6fc,afc+i]. Then, as shown by Kuijlaars and McLaughhn in 
|KuMl] . the Jacobian of the map {[oj, hi], i = 1, . . . , g} — {Tj, A^^} is nonzero at the solution, 
{[aj(t), 6i(t)], z = 1, . . . , g}, to the equations {Tj = 26jq, Nk = 0}. By the implicit function 
theorem, this implies the analyticity of [ai{t) , bi{t)] in t. □ 

20. Topological expansion 

In the paper [EM], Ercolani and McLaughlin proves topological expansion (12.151) for poly- 
nomial V of form (12.131) . with small values of the parameters tj. Their proof is based on a 
construction of the asymptotic large expansion of the parametrix for the RHP. Another 
proof of topological expansion (I2.15P is given in the paper [BI3]. Here we outline the main 
steps of the proof of [BI3]. We start with a preliminary result, which follows easily from the 
results of [DKMVZ2]. 

Proposition 20.1. Suppose V{x) is one-cut regular. Then there exists e > such that for 
all n in the interval 

n 

l-e<-<l + e, 
the recurrence coefficients admit the uniform asymptotic representation, 

7n = 7(^)+0(Ar-i), /3. = /3(^)+0(A^-i). (20.1) 
The functions 7(s), /?(s) are expressed as 

^(,) = ^(£) + K£), (20.2) 

where [a{s), b{s)] is the support of the equilibrium measure for the polynomial s^^V{x). 

The next theorem gives an asymptotic expansion for the recurrence coefficients. 

Theorem 20.2. Suppose that V{x) is a one-cut regular polynomial. Then there exists e > 
such that for all n in the interval 

l-e<^<l + e, 

the recurrence coefficients admit the following uniform asymptotic expansion as N ^ oo: 

oo 

^"~^(^)+E^'-='=/-(^). 

I.'"' o. . , (20.3) 
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where f2k{s), g2k{s), k > 1, are analytic functions on [1 — e,l + e]. 

Sketch of the proof of the theorem. The Riemann-Hilbert approach gives an asymptotic 
expansion in powers of A^~^. We want to show that the odd coefficients vanish. To prove this, 
we use induction in the number of the coefficient and the invariance of the string equations, 

7n[V'{Q)Un-l = ^ , [V'{Q)U = 0, (20.4) 

with respect to the change of variables 

H ^ l2n-j, Pj ^ /?2n-,-l, J = 0, 1, 2, . . .}. (20.5) 

This gives the cancellation of the odd coefficients, which proves the theorem. 

The main condition, under which the topological expansion is proved in [BIS], is the 
following: 

Hypothesis R. For alH > 1 the polynomial TtV{z) is one-cut regular. 

Theorem 20.3. (See [BI3].) If a polynomial V{z) satisfies Hypothesis R, then its free energy 
admits the asymptotic expansion, 

Fn - F^'''''' ~ F + Ar-2F(2) + Ar-4ir(4) + . . . ^ (20.6) 
The leading term of the asymptotic expansion is: 



1 -r 



27V)+4f(r)/3^(r)-^ 



dr, (20.7) 



where 

,ir) = ^(,) = (20.8) 

and [ci(r), 6(t)] is the support of the equilibrium measure for the polynomial V{z; r). 

To prove this theorem, we substitute asymptotic expansions (120.31) into formula ( (118.131) 
and check that the odd powers of A^~^ cancel out. See [BI3]. 

21. One-sided analyticity at a critical point 

According to the definition, see fl4.17l) - fl4.18l) . the equilibrium measure is singular in the 
following three cases: 

(1) h{c) = where c G (oj, bj), for some 1 < j < q, 

(2) h{aj) = or h{bj) = 0, for some I < j < q, 

(3) for some c ^ J, 

2 y log |c - y\dMy) - Vic) = I. (21.1) 

More complicated cases appear as a combination of these three basic ones. The cases (1) 
and (3) are generic for a one-parameter family of polynomials. The case (1) means a split of 
the support interval {aj,bj) into two intervals. A typical illustration of this case is given in 
Fig. 2. Case (3) means a birth of a new support interval at the point c. Case (2) is generic 
for a two-parameter family of polynomials. 
Introduce the following hypothesis. 

Hypothesis Sq. V{z]t), t G [0,to]) > 0, is a one-parameter family of real polynomials 
such that 
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(i) V{z; t) is g-cut regular for < t < to, 

(ii) V{z;Q) is g-cut singular and h{ai) ^ 0, h{hi) 7^ 0, i = 1, . . . , g, where U^^]^[aj,6j] is 
the support of the equilibrium measure for V{z\ 0). 

We have the following result, see [BI3]. 

Theorem 21.1. Suppose V{z]t) satisfies Hypothesis Sq. Then the end-points ai{t),bi(t) of 
the equilibrium measure for V{z;t), the density function, and the free energy are analytic, 
as functions oft, att = 0. 

The proof of this theorem is an extension of the proof of Theorem 119. 11 and it is also based 
on the work of Kuijlaars and McLaughlin |KuMlj . Theorem shows that the free energy can 
be analytically continued through the critical point, t = 0, if conditions (i) and (ii) are 
fulfilled. If h{aj) = or h{bj) = 0, then the free energy is expected to have an algebraic 
singularity at t = 0, but this problem has not been studied yet in details. As concerns the 
split of the support interval, this case was studied for a nonsymmetric quartic polynomial in 
the paper of Bleher and Eynard [BE]. To describe the result of [BE], consider the singular 
quartic polynomial, 

V^{x) = ^{x'-AciX^ + 2c2X + 8c,), T, = l + 4c?; K(0) = 0, (21.2) 
where we denote 

Cfc = cos fevre . (21.3) 

It corresponds to the critical density 

Pc{x) = ^(x - 2c,)'VI^. (21.4) 

Observe that < e < 1 is a parameter of the problem which determines the location of the 
critical point, 

- 2 < 2ci = 2cos7re < 2. (21.5) 

We include into the one-parameter family of quartic polynomials, {V{x;T), T > 0}, 
where 

V'{x; T) = ^ {x^ - 4cix2 + 2c2X + 8ci) ; V{0;T)=0. (21.6) 
Let F{T) be the free energy corresponding to V{x; T). 

Theorem 21.2. The function F{T) can be analytically continued through T = Tc both from 
T > Tc and from T < Tc. AtT = Tc, F{T) is continuous, as well as its first two derivatives, 
but the third derivative jumps. 

This corresponds to the third order phase transition. Earlier the third order phase tran- 
sition was observed in a circular ensemble of random matrices by Gross and Witten [GW]. 

22. Double Scaling Limit of the Free Energy 



Consider an even quartic critical potential, 

V{z) = \z'-z\ 



(22.1) 
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and its deformation, 

nV{z) ^ V{z- t) = ^^z'+{l--^z\ (22.2) 

Introduce the scaling, 

t = 1 + A^-2/32-2/32, _ (^22.3) 
The Tracy- Widom distribution function defined by the formula 



Ft^n{x) = exp 



(x - y)u^{y)dy 



(22.4) 



where uiy) is the Hastings-McLeod solution to the Painleve II, see (114.161) . (114.171) . 
Theorem 22.1. (See [BI3].) For every e > 0, 

FN{t) - F^^'^'' = F'^^it) - N^^ logFTw((t - l)22/3iv2/3j + 0{N-'/^+'), (22.5) 

as N ^ oo and \ {t - 1)N'^^^\ < C, where 

F'^'^'it) = F{t) + N-^F^^\t) (22.6) 
is the sum of the first two terms of the topological expansion. 

Lecture 5. Random matrix model with external source 

23. Random matrix model with external source and multiple orthogonal 

polynomials 

We consider the Hermitian random matrix ensemble with an external source, 

d^^(M) = ^e-^^-^^^^^^-^^^^rfM, (23.1) 



where 



J ^-nTrivm-AM)^^^ ^23.2) 



and A is a fixed Hermitian matrix. Without loss of generality we may assume that A is 
diagonal, 

A = diag(ai, . . . , a„), Oi < • • • < (23.3) 
Define the monic polynomial 

Pr,{z) = J det{z - M) dHn{M). (23.4) 
Proposition 23.1. (a) There is a constant Zn such that 

-. P n n 

Pn{z) = ^ \\{z- A,) n e-(^(^^)-'^^^^) A(A)rfA, (23.5) 

where 

A(A) = n(^^ - (23.6) 

and dX = d\id\2 ■ ■ ■ dXn- 
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(b) Let 



rrtjk 



Then we have the determinantal formula 

mio mil 



1 



1 z 



X. 

min 



(23.7) 



(23.8) 



1, . . . ,ri, 



(c) For j 

/oo 
P„(x)e-(^(^)-''^^)rfx = 0, (23.9) 
■oo 

and these equations uniquely determine the monic polynomial Pn- 

Proposition 123. ll can be extended to the case of multiple Ctj S clS follows. 

Proposition 23.2. Suppose A has distinct eigenvalues ai, i = l,...,p, with respective 
multiplicities Ui, so that ni + ■ ■ ■ + Up = n. Let n*^*) = ni + ■ ■ ■ + and n^^^ = 0. Define 



Wj{x) = xd,~l^^{V(x)-a^x)^ J = 1, . . . ,n, 



(23.10) 



where i = ij is such that n*^* < j < rS^^ and dj = j — rS"^ ■ Then the following hold. 



(a) There is a constant Zn > such that 



PJz) 



-. » n n 

T / ^(^-^^■)^^^•(^^■)^(^)^^• 
^n J =1 =1 



(b) Let 



rrijk 



x^Wj{x)dx. 



(23.11) 



(23.12) 



Then we have the determinantal formula 

mio mil 

PJz) ^ 



1 



nino m„i 
1 z 



min 



nir. 



(23.13) 



fc) For 2 = 1, 



oo 



P„(x)a;^e-(^(")-"'")da; = 0, j = 0, 



- 1, 



(23.14) 



and these equations uniquely determine the monic polynomial Pn 



The relations (123. 14p can be viewed as multiple orthogonality conditions for the polynomial 
Pn- There are p weights e~(^(^')~''j^\ j = 1, . . . ,p, and for each weight there are a number 
of orthogonality conditions, so that the total number of them is n. This point of view is 
especially useful in the case when A has only a small number of distinct eigenvalues. 
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23.1. Determinantal Formula for Eigenvalue Correlation Functions. P. Zinn- Justin 

proved in [ZJl], [ZJ2] a determinantal formula for the eigenvalue correlation functions of 
the random matrix model with external source. We relate the determinantal kernel to the 
multiple orthogonal polynomials. 
Let E„ be the collection of functions 

:= {x^e"'^ I i = 1, . . . ,p, J = 0, . . . - 1}. (23.15) 

We start with a lemma. 

Lemma 23.3. There exists a unique function Qn-i in the linear span ofEn such that 

/oo 
x^Qn-i{x)e-^^''^dx = 0, (23.16) 
-oo 

j — 0, . . . ,n — 2, and 

coo 

x''-^Qn-i{x)e-^'^''Ux = 1. (23.17) 

Consider Pq, . . . , P„, a sequence of multiple orthogonal polynomials such that deg — k, 
with an increasing sequence of the multiplicity vectors, so that ki < li, i = 1, . . . ,p, when 
k < I. Consider the biorthogonal system of functions, {Qk{x), k — 0, . . . ,n — 1}, 



f 

J —c 



Pj{x)Qk{x)e-^^''^dx = Sjk, (23.18) 
ior j, k — 0, . . . ,n — 1. Define the kernel 

n-l 

K^{x,y) = e-^(^(-)+^(^)) ^P,(x)g,(y). (23.19) 

k=0 

Theorem 23.4. The m-point correlation function of eigenvalues has the determinantal form 

Rmi^l, ■ ■ ^m) — det(i^^„(Aj, Xk))l<j,k<m (23.20) 

23.2. Christoffel-Darboux formula. We will assume that there are only two distinct 
eigenvalues, ai = a and 02 = —a, with multiplicities ni and n2, respectively. We rede- 
note Pn by Pni,n2- Set 

/oo 

P„„„,(x)x"^«;,(x)dx, (23.21) 
-00 

j = 1,2, which are non-zero numbers. 

Theorem 23.5. With the notation introduced above, 

_ ^)e5(^W+n.))i^^(^,^) = P,^^,^{x)Qr,,,nM - -j^Pn,-l,n,{x)Qr.,+l,nM 

^^^^ m-l,n2 ^23.22) 

~ Pni,n2—li-^)Qni,n2+liy) 



ni,n2— 1 
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23.3. Riemann-Hilbert problem. The Rieman-Hilbert problem is to find F : C \ 
C^^^ such that 

• y is analytic on C \ M, 

• for a; G M, we have 

(1 wi{x) W2{xy 
1 I (23.23) 
1 



half-plane, 

a.s z cxD, we have 



) as 










-ni 











Y{z) = [I + { - ] ] I ^-"1 I (23.24) 



where / denotes the 3x3 identity matrix. 
Proposition 23.6. There exists a unique solution to the RH problem, 

I P„,,„2 C(P„,,„2^l) C(P„,,„2W2) \ 

Y = CiP„,_i,„2 CiC{Pn^^l,n2Wl) CiC(P„^_i „2^2) (23.25) 
\C2Pni,n2-l C2C(P„i,„2-l«^l) C2C {Pm,n2~lW2) ) 

with constants 

ci = -21X1 (h'^S~i,n^ ' C2 = -27X1 , (23.26) 

and where C f denotes the Cauchy transform of f , i.e., 

Cf{z) = ^ [ I^ds. (23.27) 
2m J^s- z 

The Christoffel-Darboux formula, (123. 22p . can be expressed in terms of the solution of RH 
Problem: 

KJ., y) = e-j(n-'^v(.'>-°"^''''^"^''-""+-"'y"'^'^'<-)'^' . (23.28) 

2Txi{x — y) 

23.4. Recurrence and differential equations. The recurrence and differential equations 
are nicely formulated in terms of the function 

/I \ (w{z) \ 

m{z) = q 1 Y{z) e"^'*^ , (23.29) 

\0 V e^'^y 

where 

w{z) = e-^^(^). (23.30) 
The function \1' solves the following RH problem: 
• \E' is analytic on C \ M, 
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• for x e R, 



as z — > oo, 



A 1 l\ 

= 1 . 

\0 1/ 



(23.31) 









X ( cr^^-"ie-^«^ 
\ C2^-z""^e^' 

The recurrence equation for ^ has the form: 

*ni+l,n2(^) = t4i,n2(^)^ni,n2(^)) 

where 

V 1 

and 



(23.32) 








'ni,n2 





•ni,n2 



(1) 



_ ii'ni,n2 /p. , 
ni— 1,712 



'"1,12 ~ 772) ^ 
ni,n2 — 1 



(2) 



'ni+l,n2— 1 

•ni,n2 — (2) 

ni,n2 — 1 



(23.33) 
(23.34) 

0. (23.35) 



ni— 1,712 ni,n2 — 1 

Respectively, the recurrence equations for the multiple orthogonal polynomials are 

-fni+l,n2 (-^) — (-^ ^ni,n2)-fni,n2('^) C^i,n2-fni — l,n2 ('^) ^ni,n2-fni,n2— 1(-^) ) 

finrl 

ni+l,n2— 1 ('^) — -^ni,n2(-^) ~l~ C^i,n2-fni,n2— 1 (-^) ■ 

for ^ is 



The differential equation for ^ 



where 

^ni,n2('^) 



<i,n2(^) =^^ni,n2(^)*ni,n2(^), 



(23.36) 
(23.37) 
(23.38) 



■+^^ + ... 



'r{z) o\ / ^(1) 

7+^™ 
0/ \ ^ 



-1 



+ 



pol 



^0 
+ I -a 

lO a 



(23.39) 



where [/(2;)]poi means the polynomial part of f{z) at infinity. 

2 

For the Gaussian model, V(x) — , the recurrence eauatic 



ni+l,n2 



ni 


_!i2 


n 


n 











-2a 



1 U U I 



(23.40) 
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where n — rii + n2, and the differential equation reads 

ni n2\ 

- 1 -a (23.41) 

- 1 a J 

In what follows, we will restrict ourselves to the case when n is even and 

n 

ni = 712 = - , (23.42) 

so that 

A — diag(— a, . . . , —a, a, . . . ,a). (23.43) 



24. Gaussian matrix model with external source and non-intersecting 

Brownian bridges 



1.5r 



■ 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Figure 10. Non-intersecting Brownian paths that start at one point and 
end at two points. At any intermediate time the positions of the paths are 
distributed as the eigenvalues of a Gaussian random matrix ensemble with 
external source. As their number increases the paths fill out a region whose 
boundary has a cusp. 

Consider n independent Brownian motions (Brownian bridges) Xj{t), j = 1, . . . , n. on the 
line, starting at the origin at time t = 0, half ending at a; = 1 and half at a; = —1 at time 
t — 1, and conditioned not to intersect for t e (0, 1). Then at any time t e (0, 1) the positions 
of n non-intersecting Brownian bridges are distributed as the scaled eigenvalues. 



of a Gaussian random matrix with the external source 
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Figure 1 1 . The structure of the Riemann surface for the equation (125.11) for 
the values a > 1 (left), a = 1 (middle) and a < 1 (right). In all cases the 
eigenvalues of M accumulate on the interval(s) of the first sheet with a density 
given by (125:21) . 



FigHO] gives an illustration of the non-intersecting Brownian bridges. See also the paper 
|TW4] of Tracy and Widom on nonintersecting Brownian excursions. 

In the Gaussian model the value a = 1 is critical, and we will discuss its large n asymptotics 
in the three cases: 

(1) a > 1, two cuts, 

(2) a < 1, one cut, 

(3) a = 1, double scaling limit. 

In the picture of the non-intersecting Brownian bridges this transforms to a critical time 
t = |, and there are two cuts for t > i, one cut for t < and the double scaling limit 
appears in a scaled neighborhood of t = |. 

25. Gaussian model with external source. Main results 

First we describe the limiting mean density of eigenvalues. The limiting mean density 
follows from earlier work of Pastur |Pas] . It is based on an analysis of the equation (Pastur 
equation) 

e-ze + il-a')^ + ah = 0, (25.1) 
which yields an algebraic function C,{z) defined on a three-sheeted Riemann surface. The 
restrictions of ^{z) to the three sheets are denoted by ^j{z), j = 1,2,3. There are four real 
branch points if a > 1 which determine two real intervals. The two intervals come together 
for a = 1, and for < a < 1, there are two real branch points, and two purely imaginary 
branch points. Fig{TT] depicts the structure of the Riemann surface ^{z) for a > 1, a = 1, 
and a < 1. 

In all cases we have that the limiting mean eigenvalue density p{x) = p(x; a) is given by 

p(x; a) = -Im^i+(a;), x eR, (25.2) 

where ^i+{x) denotes the limiting value of ^i( with Im^; > 0. For a = 1 the 

limiting mean eigenvalue density vanishes at x = and p{x; a) ~ as x — 0. 

We note that this behavior at the closing (or opening) of a gap is markedly different from 
the behavior that occurs in the usual unitary random matrix ensembles Z~^e~"''^^^^^^dM 
where a closing of the gap in the spectrum typically leads to a limiting mean eigenvalue 
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density p that satisfies p(x) ~ (x — x*Y as x — x* if tlie gap closes at x = x*. In tliat case 
tlie local eigenvalue correlations can be described in terms of ■^/'-functions associated with the 
Painleve II equation, see above and [BI21 ICKj . The phase transition for the model under 
consideration is different, and it cannot be realized in a unitary random matrix ensemble. 

Theorem 25.1. The limiting mean density of eigenvalues 

p{x) = lim —Kn{x,x) (25.3) 



exists for every a > 0. It satisfies 



p(x) = - |Im^(x)| , (25.4) 



TT 

where ^ = ^(x) is a solution of the cubic equation, 

- x^2 - (a^ - 1)^ + xa^ = 0. (25.5) 

The support of p consists of those x G M for which (125.51) has a non-real solution. 

(a) For < a < 1, the support of p consists of one interval [—zi, zi], and p is real analytic 
and positive on {—zi,Zi), and it vanishes like a square root at the edge points ±zi, 
i.e., there exists a constant pi > such that 

p(x) = — |x =1= 2;i|"^/^(l + o(l)) as X ±zi, X & {—zi, zi). (25.6) 

TT 

(b) For a = 1, the support of p consists of one interval [—Zi,Zi], and p is real analytic 
and positive on {—zi, 0) U (0, zi), it vanishes like a square root at the edge points ±zi, 
and it vanishes like a third root at 0, i.e., there exists a constant c > such that 

p(x) =c|x|^/^(l + o(l)), asx^O. (25.7) 

(c) For a > 1, the support of p consists of two disjoint intervals [—Zi, —Z2] U [z2, Zi] with 
< Z2 < Zi, p is real analytic and positive on {—Zi, —z^) U {z2,zi), and it vanishes 
like a square root at the edge points ±Zi, ±^2- 

To describe the universality of local eigenvalue correlations in the large n limit, we use a 
rescaled version of the kernel Kn- 

knix, y) = e"('^(^)-'^(^»ir„(x, y) (25.8) 

for some function h. The rescaling is allowed because it does not affect the correlation 
functions Rm, which are expressed as determinants of the kernel. The function h has the 
following form on {—Zi,—Z2)l->{z2,Zi): 

1 r 

h{x) = — x^ + Re / ^i+{s)ds, (25.9) 
4 Jzi 

where ^1 is a solution of the Pastur equation. The local eigenvalue correlations in the bulk of 
the spectrum in the large n limit are described by the sine kernel. The bulk of the spectrum 
is the open interval {—zi, zi) for a < 1, and the union of the two open intervals, (— zi, —Z2) 
and {z2, zi), for a > 1 [z2 = for a = 1). 
We have the following result: 



LECTURES ON RANDOM MATRIX MODELS 



61 



\ / 

Imz 


\ 


E 









Figure 12. The contour E that appears in the definition of q{y). 



Theorem 25.2. For every Xq in the hulk of the spectrum we have that 

1 - / u V \ sin vrf-u — t>l 

hm — - — -Kn xo H -. — -, xo 



np[Xo) \ np[Xo) np[XQ) J t:[u — v) 

At the edge of the spectrum the local eigenvalue correlations are described in the large n 
limit by the Airy kernel: 

Theorem 25.3. For every u,v ^ we have 

1 / u V \ Ai(u)Ai'(v) - Ai'(u)Ai(v) 



n-*oo (p]^n)2/^ \ [piuy/^^ {piny/-^ J u — v 

A similar limit holds near the edge point —Zi and also near the edge points ±Z2 if a > 1. 

As is usual in a critical case, there is a family of limiting kernels that arise when a changes 
with n and a— i>lasn^ooina critical way. These kernels are constructed out of 
Pearcey integrals and therefore they are called Pearcey kernels. The Pearcey kernels were 
first described by Brezin and Hikami |BH4[ IBH5j . A detailed proof of the following result 
was recently given by Tracy and Widom [TW3j . 

Theorem 25.4. We have for every fixed 6 G M, 

lim -^KJ-^,-l--l + -^]=K'^-^^{x,y;b) (25.10) 



™n3/4 " ^^3/4' ^3/4' 2v/^ 

where K'^^'^^ is the Pearcey kernel 

^cusp,^ y. ^ V{x)q"{y)-p\x)q'{y) +p"{x)q{y) - bp{x)q{y) 

x-y 

with 

p(x) = — / e~i^*-t^'+'^"rfs and q{y) = — / e^^'+l^'+^'^'tit. (25.12) 

The contour S consists of the four rays argy = ±7r/4, ±37r/4, with the orientation shown in 
Fig^ 
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The functions (125. 12p are called Pearcey integrals [Peaj . They are solutions of the third 
order differential equations p'"(x) = xp{x) + bp' (x) and q"'{y) = —y(l{y) + bq'{y), respectively. 

Theorem 125.41 implies that local eigenvalue statistics of eigenvalues near are expressed in 
terms of the Pearcey kernel. For example we have the following corollary of Theorem 125. 4[ 

Corollary 25.5. The probability that a matrix of the ensemble (123.11) . (123.431) . with a = 
1 + 6n~^/^/2 has no eigenvalues in the interval [cn~^^^, dn~^^^] converges, as n ^ oo, to the 
Fredholm determinant of the integral operator with kernel K'^'^'^p{x, y; b) acting on L'^{c, d). 

Similar expressions hold for the probability to have one, two, three, . . . , eigenvalues in an 
0{n~^^^) neighborhood of a; = 0. 

Tracy and Widom |TW3j and Adler and van Moerbeke [AvM3j gave differential equations 
for the gap probabilities associated with the Pearcey kernel and with the more general 
Pearcey process which arises from considering the non-intersecting Brownian motion model 
at several times near the critical time. See also |ORj where the Pearcey process appears in 
a combinatorial model on random partitions. 

Brezin and Hikami and also Tracy and Widom used a double integral representation for 
the kernel in order to establish Theorem 125.41 We will describe the approach of [BK4j, based 
on the Deift/Zhou steepest descent method for the Riemann-Hilbert problem for multiple 
Hermite polynomials. This method is less direct than the steepest descent method for 
integrals. However, an approach based on the Riemann-Hilbert problem may be applicable 
to more general situations, where an integral representation is not available. This is the case, 
for example, for the general (non-Gaussian) unitary random matrix ensemble with external 
source, (123. ip . with a general potential V. 

The proof of the theorems above is based on the construction of a parametrix of the RHP, 
and we will describe this construction for the cases a > 1, a < 1, and a = 1. 

26. Construction of a parametrix in the case a > 1 

Consider the Riemann surface given by equation (125.11) for a > 1, see the left surface on 
FigJTTl There are three roots to this equation, which behave at infinity as 

^iW = ^-; + o(^)' =±^ + ^ + ^(^)- ^26.1) 

We need the integrals of the ^-functions, 

Xkiz) = J ^k{s)ds, k = 1,2,3, (26.2) 

which we take so that Ai and A2 are analytic on C \ (—00, zi] and A3 is analytic on C \ 
{—00, —Z2]. Then, as 2; — 00, 

Xiiz) = ^-\nz + h + o(^-^y X2,siz) =±az + Unz + l2,3 + o(^^y (26.3) 

where li, I2, I3 are some constants, which we choose as follows. We choose li and I2 such 
that 

Ai(^i) = A2(^i) = 0, (26.4) 

and then ^3 such that 

A3(-^2) = Ai+(-^2) = Ai_(-^2) - (26.5) 
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First Transformation of the RH Problem. Using the functions Xj and the constants Ij, 
j = 1, 2, 3, we define 



T{z) = diag (e-"'^e-"'^e-"'«) F(^)diag (^t 



n(Xi{z)-^z^) n{\2iz)-az) n(X3{z)+az) 



(26.6) 



Then T+(x) = T^{x)jT{x), x eM., where for x E [z2, Zi] 



JT 








gn(Ai-A2)- 





3n(A3-Ai_)^ 



1 



(26.7) 



and for x G [— ^i, —Z2] 



jT 



. gn{Ai-A3)+ gn(A2+-Ai_) ^ 

1 











(26.8) 



Second transformation of the RH Problem: opening of lenses. The lens structure is shown 
on FigJT3l Set in the right lens, 





Figure 13. The lenses for a > 1. 



S{z) 



1 

T(^z) \ — e"('^i(^)"'^2{z)) I _Qri{Xa{zy\2{z)) 

1 

in the upper lens region, 

/ 1 

\ 1 

in the lower lens region. 



(26.9) 
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and, respectively, in the left lens. 



S{z) = { 



(I 0> 

T{z) 10 

\_^n{\\(z)-\z(z)) _gn(A2(2)-A3(z)) 

in the upper lens region, 

/I 0' 

T{z) 10 

\ gn(Ai(2)-A3(z))) _gn(A2(z)-A3(2)) ^ 

in the lower lens region. 



(26.10) 



Then 



and 



S+{x) = S-{x)3s{x)\ js{x) 



1 0> 

-10 
1. 



X e [Z2,zi\, 



1' 

S+{x) ^ S4x)js{x); js{x) =10 1 0|, [-^1,-^2]. 



(26.11) 



(26.12) 



^-1 0^ 

In addition, S has jumps on the boundary of the lenses, which are exponentially small away 
of the points ±2:1 2. The RH problem for S is approximated by the model RH problem. 
Model RH problem 

• M is analytic on C \ {[—zi, —Z2] U [z2, zi]), 



M+{x) = M_{x)js{x), X e {-zi, -Z2) U (Z2, zi). 



• as 2; — > 00, 



M{z) ^I + O 



(26.13) 
(26.14) 



where the jump matrix is 



Solution to the model RH problem has the form: 







1 













< 


















\ 





1 














X e {Z2,zi) 



(26.15) 



M{z) = A{z)B{z)C{z), (26.16) 

where 

/e?(^)-a' e|(^)-a2 
A{z) = diag ( 1, ) , B(^) = I ei(^) + a 6(;^) + a U^) + a | (26.17) 

a 
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and 



C{z) = diag 



1 



(26.18) 



yVoiM^y VoiM^y VoiM^), 

where 

Q(z) =/-(! + 2a^)z^ + (a^ - l)al (26.19) 

Parametrix at edge points. We consider small disks D{±Zj,r) with radius r > and 
centered at the edge points, and look for a local parametrix P defined on the union of the 
four disks such that 

• P is analytic on D{±Zj, r) \ (M U T), 



P+{z) = P_{z)js{z), ze{RUT)nD{±Zj,r), 



as n — > oo, 



(26.20) 



P{z) = ^/ + O (^^^ M{z) uniformly for z G dD{±Zj,r). (26.21) 



II 




III 



IV 



Figure 14. Partition of a neighborhood of the edge point. 
We consider here the edge point zi in detail. We note that as z ^ zi, 
X,{z) = q{z - z^) + ^{z - z,f/' + 0{z - z,f 



so that 



\,{z) = q{z - z,) - ^-^{z - z^fl-" + 0{z - z^)^ 
X,{z) - \,{z) = ^{z - z,fl' + 0{z - z,fl' 



as 2; — > ^1. Then it follows that 



-{\,{z) - \,{z)) 



2/3 



(26.22) 



(26.23) 



(26.24) 



is analytic at Zi, real- valued on the real axis near zi and (5'{zi) = pl^^ > 0. So /3 is a 
conformal map from D{zi, r) to a convex neighborhood of the origin, if r is sufficiently small 
(which we assume to be the case). We take T near zi such that 

P{T n D{zi,r)) C {z I arg(^) = ±27r/3}. 



66 



PAVEL M. BLEHER 



Then F and M divide the disk D{zi,r) into four regions numbered I, II, III, and IV, such 
that < aig /3{z) < 2tt/3, 2tc/3 < &Tgf3{z) < tt, -tt < &rg /3{z) < -2tt/3, and -27r/3 < 
arg/3(z) < for 2; in regions I, II, III, and IV, respectively, see FigJT^ 
Recall that the jumps js near zi are given as 



Js 



3s 



Js 




on [Zi — r,Zi] 



1 

gn(Ai-A2) ]^ gn(A3-A2) on the uppcr bouudary of the leus in D (zi , r ) 











1 

gn(Ai-A2) ]^ _gn{A3-A2) on the lowcr bouudary of the Icus in r) 



(26.25) 











]s = 
We write 



'1 g"{A2-Ai) gn(A3-Al)^ 

1 
1 



on (zi, zi + r] 



P 






3n(A3-A2) 
1 



in regions I and IV 
in regions II and III. 



(26.26) 



Then the jumps for P are P+ = P-jp where 



Jp 



Jp 



Jp 



Jp 




on zi -r,zi, 



1 0^ 

gra(Ai-A2) 1 on the upper side of the lens in D(zi,r) 







1^ 
0^ 



(26.27) 



^n{Xi A2) ^ g I j^Yie lower side of the lens in D{zi,r) 
1, 



'1 e"*^'^^"'^^-* 0^ 
1 
,0 



on [Zi, zi + r . 



We also need the matching condition 



n 



P[z) = / + O - M{z) uniformly for z G dD{zi,r) 



(26.28) 



The RH problem for P is essentially a 2 x 2 problem, since the jumps (126.271) are non-trivial 
only in the upper 2x2 block. A solution can be constructed in a standard way out of Airy 
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functions. The Airy function Ai{z) solves the equation y" = zy and for any e > 0, in the 
sector 7T + e < aigz < tt — e, it has the asymptotics as z oo, 

1 



M{z) 



_2 3/2 

-e 3^ 



(1 + 0(^-3/2)) 



(26.29) 



2v/^zV4- 

The functions Ai(ujz), Ai(uj'^z), where lu = e^, also solve the equation y" = zy, and we 
have the linear relation, 

Ai(;z) + ujAi{ujz) + uj^Ai{uj^z) = 0. (26.30) 

Write 

yo{z) = Ai{z), yi{z) = ujAi{ujz), y2{z) = uj^Ai{uj^z), (26.31) 
and we use these functions to define 



'y,{z) -y<,{z) 0^ 
y'^{z) -y'^iz) 
1^ 

-y^{z) -y2{z) 0' 
-y[{z) -y'^{z) 
1 

M^) yi{z) 0^ 

1, 



for < argz < 27r/3, 

for 27r/3 < axgz < n, 
for —71 < aTgz < — 27r/3, 



(26.32) 



/yoiz) y,{z) 0' 

y'Q^z) y[{z) I , for -27r/3 < argz < 0. 
\ 



Then 



P{z) = E„(z)<l>(n2/3/5(;z))diag (^e'^M^'^-)->^2{z)) ^ ^~'^n{x,{z)-x,{z)) ^-^-^ 



(26.33) 



where En is an analytic prefactor that takes care of the matching condition (126.281) . Explic- 
itly, En is given by 



En = v^M 




nV6/5i/4 0^ 

n-i/6/5-1/4 I . 
1, 



(26.34) 



A similar construction works for a parametrix P around the other edge points. 
Third transformation. In the third and final transformation we put 

R{z) = S{z)M{z)-'^ for z outside the disks D{±Zj, r), j = 1, 2 
R{z) = S{z)P{z)^'^ for z inside the disks. 



(26.35) 



Then R is analytic on C \ Tn, where consists of the four circles dD{±Zj,r), j = 1,2, 
the parts of T outside the four disks, and the real intervals {—oo, —zi — r), {—Z2 + r, Z2 — r), 
{zi + r, oo), see Fig|T3 There are jump relations 



R+ = R^jR 



(26.36) 
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Figure 15. The contour Tji for R. 



where 



jji = MP ^ on the circles, oriented counterclockwise 
jji = MisM~^ on the remaining parts of Vji. 



(26.37) 



We have that = I + 0{l/n) uniformly on the circles, and jji = I + 0{e '^") for some c > 
as n ^ oo, uniformly on the remaining parts of F/j. So we can conclude 



n 



as n — >■ cxo, uniformly on F/j. 



As 2; — > oo, we have 



R{z) = I + 0{l/z). 



(26.38) 



(26.39) 



From (126. 36p , (126. 38p , (126.391) and the fact that we can deform the contours in any desired 
direction, it follows that 



R{z) =1 + 



1 



as ri ^ oo. 



n{\z\ + 1 

uniformly for 2 G C \ Fr, see IDKMVZl IDKMVZ21 Km]. 



(26.40) 



27. Construction of a parametrix in the case a < 1 

27.1. A-functions. Consider the Riemann surface given by equation (125.11) for a < 1, see 
the right surface on FiglTTl There are three roots to this equation, which behave at infinity 
as in (126. ip . We need the integrals of the ^-functions, which we define as 

-^i(^) = / ^i{s)ds, 

J Zl 

Hz) = I Us)ds, (27.1) 

J Zl 

\^{z) = I i:i{s)ds + Xi_{-Zi), 

The path of integration for A3 lies in C\ ((—00, 0] U [—iz2., iz2\), and it starts at the point —Zi 
on the upper side of the cut. All three A-functions are defined on their respective sheets of 
the Riemann surface with an additional cut along the negative real axis. Thus Ai, A2, A3 are 
defined and analytic on C\(— 00, zi], C\((— 00, zi]U[—iz2, iz2]) , and C\((— 00, Q]U[—iz2,iz2]), 
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(27.3) 



respectively. Their behavior at infinity is 

Hz) = U^-\ogz + h + 0{llz) 

\2{z)^az + \\ogz + £2 + 0{l/z) (27.2) 

\j,{z) ^-az+^ log z + i3 + 0(1/ z) 

for certain constants £j, j — 1,2,3. The Aj's satisfy the following jump relations 

Ai=F = A2± on (0,zi), 

Ai_ = A3+ on(-2;i,0), 

Ai+ = A3_-7ri on(-zi,0), 

= ^3± on (0,2^2), 

A2:f = ^3± - on {-iz2, 0), 

Ai+ = Ai_ — 27ri on (—00, —Zi), 

A2+ = A2-+7ri on(-cx),0), 

A3+ = A3_ + 7ri on(-oo, -^1), 

where the segment {—iz2,iz2) is oriented upwards. 

27.2. First transformation Y ^U. We define for z e C \ (M U [-iz2,iz2]), 

U{z) = diag(e-"^\ e""^^ e-"^3)F(z)diag(e"(^i(^)-5^'\ e^^^^W-a^)^ ^n{X;iz)+az)^_ ^^27.4) 

This coincides with the first transformation for a > 1. Then U solves the following RH 
problem. 

• : C \ (M U [-iz2, iz2]) C^""^ is analytic. 

• U satisfies the jumps 

(gn(Ai+-Ai_) gn(A2+-Ai_) gn(A3+-Ai_)\ 
^n{X2+-X2-) I onR, (27.5) 

gn(A3+-A3-) J 

and 

/I \ 

= e'*(^2+-^2-) on[-iz2,iz2]. (27.6) 

\0 f>n{X3+-X3-) J 

• U{z) = I + 0{l/z) as z^ 00. 

27.3. Second transformation U ^ T: global opening of a lens on [—iz2,iz2\. The 
second transformation is the opening of a lens on the interval [—iz2,iz2\. We consider a 
contour E, which goes first from {—iz2) to iz2 around the point Zi, and then from iz2 to 
{—iz2) around the point —zi, see Fig. 16, and such that for 2; e E, 

± (ReA2(;2) -ReA3(z)) > 0, ±Rez>Q. (27.7) 

Observe that inside the curvilinear quadrilateral marked by a solid line on Fig. 16, ± (Re A2(-2) — Re A3(^)) < 
0, hence the contour E has to stay outside of this quadrilateral. We set T = 6'^ outside E, 
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Re z 



Figure 16. Contour E which is such that ReA2 < Re A3 on the part of S in 
the left half-plane and Re A2 > Re A3 on the part of E in the right half-plane. 

and inside E we set 

/I 0\ 

T = [/ 1 for Rez < inside E, 

VO -e"(^2-A3) I 

(27.8) 

'10 0^ ^ ' 

T \ 1 -e"(^:^-^2) I for Re^ > inside E. 
^0 1 

27.4. Third transformation T S: opening of a lens on [—zi, zi]. We open a lens on 
[— -^1, -^i], inside of E, see Fig. 17, and we define S = T outside of the lens and 

1 0\ 

S = T \ ^ ^ \ upper part of the lens in left half-plane, 

gn(Ai-A3) g 

1 0' 

S — T I 10 1 in lower part of the lens in left half-plane. 



gn(Ai-A3) g ^ 

1 0\ 

S — T \ — e"('^i~'^2) 10 in upper part of the lens in right half- plane, 

1/ 

1 0\ 

S — T \ e"('^i~'^2) 1 I in lower part of the lens in right half- plane. 
1/ 

Then S satisfies the following RH problem: 



(27.9) 
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Figure 17. Opening of a lens around [—2:1, zi]. 



• S is analytic outside the real line, the vertical segment [—iz2,iz2], the curve E, and 
the upper and lower lips of the lens around [—zi, zi\. 

• S satisfies the following jumps on the real line 



S-L. S- 



— S- 



S-L. S— 



^-L S- 



= s. 



'-^ gn(A2+-Ai_) gn(A3+-Ai-) 

1 I on (-00, -X* 

^0 1 

1 I on {-x\ -zi\ 

,0 1 



1' 
1 
-1 0^ 

1 0^ 
-10 
1^ 

gn(A2-Ai) qN 





1 
,0 



'\ gn(A2-Ai) gn(A3-Ai) 



1 
.0 



on (-2:1,0) 
on (0, zi) 

on [-^i, X*) 

on [a;*, 00). 




1 



(27.10) 
(27.11) 
(27.12) 
(27.13) 
(27.14) 
(27.15) 
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5* has the following jumps on the segment [—iz2, iz2], 

/I \ 

= S_ 1 on {~iz2, ~iy*) (27.16) 

y —1 e"'^'*'^+~'^^"^ / 

/ 1 \ 

5+ = 5_ 1 on(-zy*,0) (27.17) 

lgn(Ai-A3_) _^ gn(A3+-A3_) / 

/I \ 

5+ = 5_ 1 on{0,iy*) (27.18) 

\ _gn(Ai-A3_) _^ gn(A3+-A3_) / 

on{iy*,iz2). (27.19) 

—1 g"(A3+-A3-) / 

The jumps on S are 

/I 0\ 

5+ = 5_ 1 on {z G S I Rez < 0} (27.20) 

\o e"(^2-A3) 

/I \ 

5+ = 5_ 1 e"(^3-^2) on {z G S I Rez > 0}. (27.21) 

Vo 1 ; 

Finally, on the upper and lower lips of the lens, we find jumps 
0^ 

5*+ = S*- I I the lips of the lens in the left half-plane (27.22) 





S",^ = 5'_ en(Xi~\2) I Q\ on the lips of the lens in the right half-plane. (27.23) 



• S{z) = I + 0(1/ z) as z ^ oo. 

As n ^ oo, the jump matrices have limits. Most of the limits are the identity matrix, except 
for the jumps on {-Zi,Zi), see (127. 12p and (127.131) . and on {-iz2,iz2), see ( l27.16l) - (l7r:Tg]) . 
The limiting model RH problem can be solved explicitly. The solution is similar to the case 
a > 1, and it is given by formulas (I26.16p - (I26.19I) . with cuts of the function a/ P{z) on the 
intervals [— and [—iz2,iz2]. 

27.5. Local parametrix at the branch points for a < 1. Near the branch points the 
model solution M will not be a good approximation to S. We need a local analysis near each 
of the branch points. In a small circle around each of the branch points, the parametrix P 
should have the same jumps as S, and on the boundary of the circle P should match with 
M in the sense that 

P{z) = M{z) (/ + 0{l/n)) (27.24) 
uniformly for z on the boundary of the circle. 
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The construction of P near the real branch points ±zi makes use of Airy functions and it 
is the same as the one given above for the case a > 1. The parametrix near the imaginary 
branch points -^1x2 is also constructed with Airy functions. We give the construction near 
iz2- There are three contours, parts of S, meeting at 1x2. left, right and vertical, see Fig.l7. 
We want an analytic P in a neigborhood of iz2 with jumps 



on left contour 



P^ 



P^ 




on right contour 



(27.25) 



In addition we need the matching condition (127.241) . 
(127.241) . the problem is a 2 x 2 problem. We define 



on vertical part. 



Except for the matching condition 



A3)(^) 



2/3 



(27.26) 



such that 



arg/(2;) = 7r/3, ioi z = iy, y > Z2. 

Then s = f{z) is a conformal map, which maps [0,i2;2] into the ray args = — and which 
maps the parts of S near iz2 in the right and left half-planes into the rays args = and 
respectively. The local parametrix has the form. 



args 



P{z) = E{z)<i> {n^^'^fiz)) 







A3) 




(27.27) 



e' 2' 

where E is analytic. The model matrix-valued function $ is defined as 




for < args < 27r/3, 



for — 27r/3 < arg s < 0, 



for 27r/3 < args < 47r/3, 



(27.28) 



where yo{s) = Ai(s), yi{s) = uAi^us), y2{s) = u'^Ai^u^s) with u = 27r/3 and Ai the standard 
Airy function. In order to achieve the matching (127.241) we define the prefactor E as 

E = ML-^ (27.29) 
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with 



■1 
2v^ \ Q 




(27.30) 



where /^/^ has a branch cut along the vertical segment [0, iz2] and it is real and positive where 
/ is real and positive. The matching condition (127.241) now follows from the asymptotics of 
the Airy function and its derivative. 

A similar construction gives the parametrix in the neighborhood of —iz2- 

27.6. Fourth transformation S ^ R. Having constructed and P, we define the final 
transformation by 



R{z) = S{z)M{z) ^ away from the branch points, 
R{z) = S{z)P{z)^'^ near the branch points. 



(27.31) 



Since jumps of S and A^ coincide on the interval {—zi, zi) and the jumps of S and P coincide 
inside the disks around the branch points, we obtain that R is analytic outside a system of 
contours as shown in Fig. 18. 




Figure 18. R has jumps on this system of contours. 

On the circles around the branch points there is a jump 

R+ = R^{I + 0{l/n)), (27.32) 
which follows from the matching condition (127.241) . On the remaining contours, the jump is 

i?+ = i?_(/ + 0(e~"")) (27.33) 
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for some c > 0. Since we also have the asymptotic condition R{z) = I + 0(1/ z) a.s z ^ oo, 
we may conclude as that 

R{z) = I + o(—^ -) asn^cx), (27.34) 



n{\z\ + 1) 



uniformly for 2; G C 



28. Double scaling limit at a = 1 
This section is based on the paper |BK4j . 

28.1. Modified Pastur equation. The analysis in the cases a > 1 and < a < 1 was 

based on the Pastur equation (125. ip . and it would be natural to use (125. ip also in the case 
a = 1. Indeed, that is what we tried to do, and we found that it works for a = 1, but in 
the double scaling regime a = 1 + with 6 7^ 0, it led to problems that we were unable to 
resolve in a satisfactory way. A crucial feature of our present approach is a modification of 
the equation (125.11) when a is close to 1, but different from 1. At x = we wish to have a 
double branch point for all values of a so that the structure of the Riemann surface is as in 
the middle figure of Figure [TT] for all a. 

For c > 0, we consider the Riemann surface for the equation 

z = (28.1) 

where is a new auxiliary variable. The Riemann surface has branch points at z* = ^^c, 
—z* and a double branch point at 0. There are three inverse functions Wk, k = 1,2, 3, that 
behave as 2; 00 as 

/ 1 

wi{z) = z h O — 

z yz'^ 

„,W=c+|l + 0(^i) (28,2) 

2z \z 

and which are defined and analytic on C \ [—z*, z*], C \ [0, z*] and C \ [—z*, 0], respectively. 
Then we define the modified ^-functions 

^k = Wk + —, for A; = 1,2, 3, (28.3) 
which we also consider on their respective Riemann sheets. In what follows we take 

and p = c^-l. (28.4) 



wsiz) = -c+ — + 0{ — 



4 

Note that a = 1 corresponds to c = 1 and p = 0. In that case the functions coincide with 
the solutions of the equation (125.11) that we used in our earlier works. From (128. ip . (128. 3p . 
and (128. 4p we obtain the modified Pastur equation 

e - ze + (1 - a')i + a'z + = 0, (28.5) 



76 PAVEL M. BLEHER 

where c is given by (128. 4p . This equation has three solutions, with the following behavior at 
infinity: 



1 / . X (28-6) 



6,3(^) = ±a + - + 0^^2 
and the cuts as in the middle figure of FigiHl At zero the functions ^fc have the asymptotics, 

-^'^^z^l^j2{z) - u^z-^l^g2{z) + - for Im^ > 0, 

UA={ , , I (28.7) 

-u''z^/'^f2{z) - uj^''z-^/^g2iz) + - for Imz < 0, 

where the functions f2{z),g2{z) are analytic at the origin and real for real z, with 

/2(0) = c^/^ + lc~'/\c' - 1), 92(0) = c~'/-\c' - 1). (28.8) 
We define then the functions as 

>^kiz) = [ Us)ds (28.9) 



where the path of integration starts at on the upper side of the cut and is fully contained 
(except for the initial point) in C \ (—00,2;*], and we define the first transformation of the 
RHP by the same formula (126. 6p as for the case a > 1. For what follows, observe that the 
A-functions have the following asymptotics at the origin: 

-^u^'^z'/'Mz) - lu'z'/-'g,{z) + ^ for Im z > 0, 

h{z)={ o . 2 (28-10) 

A,_(0) - ^oo'z'/'Mz) - -oo^'z'/'g.iz) + j for Imz < 0, 

where the function /s and gs are analytic at the origin and 

/3(0) = /2(0) = c'/' + lc-'/%c' - 1), gsiO) = 3g2{0) = ^c-"\c' - 1), (28.11) 




Figure 19. Lens structure for the double scaling limit. 

The second transformation, the opening of lenses is given by formulas (126.91) . (126. lOp . The 
lens structure is shown on Fig. 19. The model solution is defined as 

fM,{w,{z)) M,{w2{z)) M^{w,{z))\ 
M{z) = M2{w,{z)) M2{w2{z)) M2{ws{z)) (28.12) 
\Ms{w,{z)) M3{w2{z)) M,{ws{z))J 
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Mi(w) 



w 



w'^ — 



w + c 



w — c 



\/2 w\/w^ — 3c^ 



(28.13) 



The construction of a parametrix P at the edge points ±z* can be done with Airy functions 
in the same way as for a > 1. 



28.2. Parametrix at the origin. The main issue is the construction of a parametrix at 
the origin, and this is where the Pearcey integrals come in. The Pearcey differential equation 
P"'(C) = Cp(C) + ^P'(C) admits solutions of the form 



_is4_ b 2 , 



(28.14) 



for j 



0,1,2,3,4,5, where 



Fo = (-00, cx)). 

Fa = (200,0] U [0, -00), 

F4 = (-icx),0] U [0, CX)), 



Fi = (zcx),0] U [0,00), 
Fg = (-icx),0] U [0, -00), 
F5 = {—ioo, zoo) 



^28. 15) 



or any other contours that are homotopic to them as for example given in Fif 
formulas (128.151) also determine the orientation of the contours Fj. 



The 




Figure 20. The contours F^, j = 0, 1, . . . , 5, equivalent to those in (128.151) . 
that are used in the definition of the Pearcey integrals . 
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Define $ = &) in six sectors by 



-P2 Pi P5 

^ = I -P2 P'l P'5 

-p'2 Pi p'L 

'Po Pl Pi 
$ = I p'o p'l p'a 
^Po P'i P'i, 

-P3 -P5 Pi ' 

= I -P's -p'5 P'4 

-Pl -Pl Pl 

'Pa -P5 P3' 

$ = I p1 -pI pI 

,p1 -Pl Pl 

^PO P2 P3^ 
^ = I Po P'2 Pl 

,Pl P2 Pl 
^Pl P2 Pb 

^= I p'l p'2 pI 
p'l p'2' p'l 



for < arg^ < 7r/4 
for 7r/4 < argC < 37r/4 
for 37r/4 < arg^ < tt 

for — TT < argC < —Sir/ 4 
for - 37r/4 < argC < -n/A 
for — 7r/4 < argC < 



We define the local parametrix Q in the form 

Qiz) = EizMn^^'Ciz); v'l%{z))e^'^'^''^ e-^""" l\ A = diag(Ai, A2, A3 
where E is an analytic prefactor, and 

C,{z) = ({z; a) = z [f2,{z] a)]^'^ 

and 

b[z) = b{z; a) — 



(28.16) 
(28.17) 
(28.18) 
(28.19) 
(28.20) 
(28.21) 

(28.22) 
(28.23) 
(28.24) 



/3(^;a)i/2- 

The functions /s, appear in (128. lOp . in the asymptotics of the A-functions. In (128. 23p and 
(128.240 the branch of the fractional powers is chosen which is real and positive for real values 
of z near 0. The prefactor E{z) is defined as 

'n^/^ 



Eiz) 




^ie-^b{^r/^M{z)K{C{z))~^ 
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where 



f /(-'/^ \ /-iu uj^ 1' 



1 



^^-1/3 

1 





-1 1 1 

— tU^ UJ 1 , 
'uj"^ UJ V 

1 1 1 

UJ up' \ 



n-^l^ , 



for Im(^ > 0, 



for Im(^ < 0. 



(28.25) 



(28.26) 
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28.3. Final transformation. We fix 6 G M and let a = 1 + 7r7= and we define 



R{z) 



' S{z)M{z)-'^, for 2 e C \ Ss outside the disks D{Q,n-^'*) and D{±^ 

S{z)P{z)-\ ioi zeD{±'^,r)\Y.s. 
{ S{z)Q{zY\ for z e Z^(0, n-"^) \ S5. 



2 ''^J' 



(28.27) 

Then R{z) is analytic inside the disks and also across the real interval between the disks. 
Thus, R{z) is analytic outside the contour S/j shown in Figl2Tl On the contour E/j the 




Figure 21. The contour S/j. The matrix- valued function R is analytic on 
C \ S/j. The disk around has radius n~^/^ and is shrinking as n 00. The 
disks are oriented counterclockwise and the remaining parts of are oriented 
from left to right. 



function R{z) has jumps, so that R+{z) ■ 

jR{z) = I + 0{n-') 

jn{z)=I + Oin-'/') 
and there exists c > so that 



R-iz)jR{z), where 



uniformly for z =]= 
uniformly for \z 



3^3 



n 



-1/4 



jr{z) =1 + 



l + \z\^ 

Also, as z 00, we have R{z) 
uniformly for z G C \ Sr. 



uniformly for z on the remaining parts of S/j. 

1 + 0(1/2;). This implies that 

^-1/6 

R{z) =1 + 0' 



l + \z\ 



(28.28) 
(28.29) 

(28.30) 
(28.31) 



29. Concluding remarks 



The Riemann-Hilbert approach is a new powerful method for random matrix models and 
orthogonal polynomials. In this paper we reviewed the approach and some of its applications. 
Let us mention some recent developments. The RH approach to orthogonal polynomials with 
complex exponential weights is considered in the recent work of Bertola and Mo |BMj . The 
RHP for discrete orthogonal poljTiomials and its applications is developed in the monograph 
[BKMMj of Baik, Kriecherbauer, McLaughlin and Miller. Applications of random matrix 
models to the exact solution of the six-vertex model with domain wall boundary conditions 
are considered in the works of Zinn- Justin |ZJ3j and Colomo and Pronko [CPj . The RH 
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approach to the six-vertex model with domain wall boundary conditions is developed in 
the work of Bleher and Fokin |BF] . The RHP for a two matrix model is considered in 
the works of Bertola, Eynard, Harnad |BEH2] and Kuijl and McLaughlin [KuM2j . The 
universality results for the scaling limit of correlation functions for orthogonal and symplectic 
ensembles of random matrices are obtained in the works of Stojanovic [Stoj . Deift and Gioev 
[DGlJ,[DG2j, Costin, Deift and Gioev [CDGj . Deift, Gioev, Kriecherbauer, and Vanlessen 
[DGKVj . 
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